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Abstract Smcctic order on arbitrary curved substrate can be described by 
a differential form of rank one (1-form), whose geometric meaning is the 
differential of the local phase field of the density modulation. The exterior 
derivative of 1-form is the local dislocation density. Elastic deformations are 
described by superposition of exact differential forms. We use the formalism 
of differential forms to systematically classify and characterize all low energy 
smectic states on torus as well as on sphere. A two dimensional smectic 
order confined on either manifold exhibits many topologically distinct low 
energy states. Different states are not accessible from each other by local 
fluctuations. The total number of low energy states scales as the square root 
of the system area. We also address the energetics of 2D smectic on a curved 
substrate and calculate the mean field phase diagram of smectic on a thin 
torus. Finally, we discuss the motion of disclinations for spherical smectics 
as low energy excitations, and illustrate the interesting connection between 
spherical smectic and the theory of elliptic functions. 
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1 Introduction 

The term smectic order [8lll5j. usually refers to one dimensional modula- 
tion of certain physical quantity with well defined periodicity. Being only 
partially ordered, smectic systems exhibit solid-like rigidity only in one di- 
mension, together fluidity in the other dimensions. As a consequence, their 
responses to external stress and fields are extremely anisotropic and exotic. 
The physics of smectic order, including the anomalous elasticities result- 
ing from thermal [17U18ifl"6l l7] and quenched fluctuations [44043113] . as well 
as nematic-Smectic transitions [I4"ll2tjll55] . have been subjected to intensive 
studies in the past few decades. 

Recent studies of crystalline order on sphere i-e. the generalized 

Thompson problem, have revealed rich and striking interplays between de- 
fects, topology, and geometry. What is perhaps most surprising out of these 
studies is the fractionalization of +1/6 disclinations into grain boundaries of 
finite extension. The appearance of these string-like objects is purely due to 
the interaction between elastic strain and Gaussian curvature, and has no 
analogue in flat space. More generally, the connection between defects, ge- 
ometry and orders has been one of the man theme of soft condensed matter 
physics in the past two decades [57] . 

The problem of smectic order on curved surfaces have attracted sub- 
stantial attention from theorists and experimentalists recently. For example, 
Santangelo et. al. [ITJ have analyzed the energetics of smectic phase formed 
by block copolymer films coated on a curved substrate. Coating of colloidal 
particles with diblock copolymer films may also provide alternative mecha- 
nism for self-assembling, which can be further exploited for designs of novel 
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artificial atoms [36]. Kleman and Blanc |4] studied a colloidal particle im- 
mersed in a smectic liquid crystal. If the boundary conditions are such that 
the nematic order must be tangent to the boundary, the smectic layers in- 
tersect the surface vertically. The loci of intersection then naturally define 
a stack of equal-distance lines, i.e. two dimensional smectic order, on the 
boundary surface. Smectic orders are also believed to arise when a flexi- 
ble charged polymer is adsorbed onto an oppositely charged curved surface, 
such as that of a colloidal particle [Hfll^l3"Ill40lPfTj . As a consequence of 
the electrostatic attraction from the substrate, as well as the repulsion be- 
tween different chain segments, the polymer may form equal-distanced layer 
pattern on the surface. More recently, polymer vesicles with in-plane smectic 
order has been synthesized [28j . which provides an experimental realization 
of smectic order on two-dimensional flexible membranes. Last but not least, 
motivated by structures formations in biological systems, stripe patterns on 
sphere were simulated via reaction-diffusion systems [521132] . 

In the two dimensional flat space, dislocations of smectic orders, cost only 
finite energy, and hence always destroy the long range smectic order at any 
finite temperature. Nevertheless, our primary interest in this work is smectic 
orders on compact, i.e. closed, substrates that are of finite size and have no 
boundary, such as torus and sphere. Typical examples include liquid film 
with in-plane modulation on the surface of a colloidal particle at the micron 
scale, and a polymer vesicle with in-plane stripe pattern. We are therefore 
free of taking any thermodynamic limit. In sufficiently low temperature, then, 
proliferation of dislocations is strongly suppressed, and the relevant states are 
the low energy states with minimal number of defects. A question of primary 
importance is the existence of defects-free smectic states on a given substrate 
topology. If the answer is negative, what do minimally defected state look 
like? More generally, we are interested in classifying all low energy, minimally- 
defected states on topological basis, and finding their connections with the 
substrate topology. The answer to the above questions are of topological 
significance, and hence apply to smectic order both on rigid curved substrates 
and on flexible membranes. Simultaneously, we are also interested in the 
geometric as well as energetic aspects of low energy smectic states on rigid, 
curved substrates Q. The energetics of smectic order on flexible membranes 
is more complicated and will be addressed in a separate publication [56j . In 
the exploration of these topological and geometric questions, we shall find 
the modern formalism of differential forms especially helpful. 

A few generic results about smectic order on compact substrates emerge 
from our study. Firstly, smectic orders both on sphere and torus exhibit many 
topologically distinct and minimally defected low energy states. These states 
can be classified by two integer topological charges, which arc related to the 
substrate topology and to the disclination pattern. They are approximately 
degenerate to each other, with the energy differences scaling sublinearly with 
the system size. The total number of these low energy states scales as the 
square root of the system size. For a smectic on sphere, twisting of closely 



1 We note that some interesting results on the interplay between smectic order 
and extr insi c as well as intrinsic curvature have been explored by Santangelo et al 
recently [37] . 
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bound disclination pairs constitutes an interesting mode of low energy exci- 
tation that allows us to prove the topological equivalence of many seemingly 
very different low energy states. For smectic on a torus, we study the ener- 
getic couplings between smectic order and extrinsic curvature in more detail. 
We calculate the mean field phase diagram and find a continuous chiral- 
achiral transition that is driven by the competition between elastic strain 
energy and the energetic coupling between director field and the extrinsic 
curvature. Some parts of our results were already published elsewhere |55j . 

There are two complementary descriptions of smectic ordering. In the 
weakly ordered regime, a smectic order can be described as one dimensional 
modulation of the relevant physical quantity with preferred wave length do'. 



where 0(x) is the phase field of the density modulation. In the flat space 
and using Cartesian coordinate system, the ground state is characterized by 
a phase field 



with (jo = 2ir/do the wave vector, and the unit vector n the direction of 
modulation. In the ground state Eq. (J2J in the flat space, 0(x) is a globally 
well-defined field. This is no longer true in the presence of dislocations, or 
when the system is confined in a closed surface. Moreover, the density profile 
Eq. ([!]) is invariant under the transformation 0{x) — > —0(x), i.e. the unit 
vector n is defined only up to a global sign. This fact is intimately related to 
the 7r-rotation symmetry of the underlying ncmatic order, and leads to the 
existence of ±1/2 disclinations. 

In the strongly ordered regime, a purely geometric description is more 
economical, where the smectic order is characterized by a stack of mathe- 
matical, i.e. infinitely thin, surfaces, which correspond to the ridges of the 
density profile Eq. (JTJ>. In a d dimensional space, each surface, called a smec- 
tic layer, is a d — 1 dimensional submanifold. In the flat space and in the 
absence of disloation, it is possible to assign a distinct phase 2kir to each 
layer, where k is an integer. For example, in the ground state Eq. ([2]), each 
layer is a plane given by qoh ■ x = 2fc7r; the separation between neighboring 
layers is uniformly do = 2n/qo. Elastic deformations can be described by 
continuous, single valued phonon fields, defined relative to the ground state. 
Dislocations are traditionally described, again relative to ground state by 
singular, multiple- valued displacement fields that are properly quantized [5]. 
Equivalcntly, the phase field 0(x) in Eq. (fTJ) becomes multiple- valued in the 
presence of dislocations. 

On a compact substrate, the phase field 0(x) may be multiple- valued even 
in the absence of defects. Careful analysis of this strange possibility leads to 
characterization of all low energy smectic states with minimal number of 
defects on closed surfaces. As a simple example, let us consider a defects-free 
smectic state on a cylinder Q. As shown in Fig. [TJ we locally identify a layer 

2 A cylinder is not a compact manifold such as torus or sphere. It however serves 
as the simplest example that exhibit many topologically distinct low energy smectic 
states. 



p(x) = po + Pi COS0(x), 



(1) 



0(x) = qo h ■ x, 



(2) 
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Fig. 1 (Color online) We label different layers by integer n = 0,1,2,... on a 
cylinder. As the thick layer (black) wraps around the cylinder once and comes 
back, however it becomes the layer 3. This integer 3 behaves as a topological charge 
and classify all low energy smectic states on a cylinder with no local defects. This 
global topological charge is stable against all local fluctuations. 




A B 

Fig. 2 How to change the global dislocation charge k: A) Starting from a state 
with k = 2, cut the whole cylinder open along the cylinder axis. B) Shear the 
cylinder and reconnect all broken layers in a different way. The final state is k = 1. 
Equivalently one can also generate a pair of dislocations (layer terminations), by 
cutting one layer, and let them propagate in opposite direction through the whole 
cylinder. Either way results in a global rearrangement of the smectic layer pattern. 



and label it with integer 0. We shall also chose a positive direction and label 
all neighboring layers successively as — 1, — 2, . . . and 1,2,.... As shown in 
Fig. [U as the 0-th layer wraps around the cylinder once along the right hand 
direction, it comes back as the 3-rd layer. The phase field is therefore well 
defined only up to 6-7T. In a more general case, clearly, the phase uncertainty 
is 2fc7r, where k is an arbitrary integer. The phase field 0{x) is single- valued 
oly if k = 0. It is easy to see that the integer k is also the minimal number 
of smectic layers that one has to cut in order to cut the cylinder into two 
halves in the direction perpendicular to the axis. 

Note that two states with different k's are not continuously accessible 
from each other. For such a purpose, a global and discontinuous rearrange- 
ment is needed. For example, one may cut all layers along the direction of 
cylinder axis (Fig. [2lA_) and reconnect them in a different way (Fig. [2)3) , or, 
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equivalently, one may generate a pair of dislocations (by cutting one layer) 
and move them along opposite direction out of the system. In each process, 
one has to cut infinite number of smectic layers. Therefore on a cylinder there 
exist an infinite number of defects-free smectic states that are topologically 
distinct from each other, each characterized by an integer k. 

As will be discussed in detail in Sec. [H the total free energy of smectic 
order on curved substrate consists of the energy cost for smectic layer com- 
pression/dilation and the layer bending energy cost. It is clear that for all 
states the layer spacing can be continuously tuned to be do, the energeti- 
cally preferred value. Therefore the layer compression/dilation energy cost 
vanishes for every k. The bending energy cost, on the other hand, is different 
for each k, and gencrically depends on several clastic moduli. Nevertheless, 
its order of magnitude can be estimated as /bending = -fbendingA4 ~ K/ R 2 , 
where L and R arc the height and radius of the cylinder respectively. The free 
energy difference per unit area of all states with different integer k therefore 
scales inversely with R 2 : all these topologically distinct smectic states on a 
cylinder become approximately degenerate as the system size becomes large 
comparing to the layer spacing. 

An integer that classifies degenerate but topologically distinct ground 
states is usually called a topological quantum number [511154] . We shall see 
that the integer k in the present case also shares essential similarities with the 
Burgers' number of a dislocation, but originates from the global (i.e. topo- 
logical) structure of the underlying substrate. It therefore may also be called 
a global dislocation charge. As we shall show in great detail in this work, the 
existence of many topological distinct low energy states is a generic property 
of smectic on compact substrates. In contrast, all defects free smectic states 
in flat space are topologically identical to each other. 

It is important to note that crystalline orders on cylinder has been inten- 
sively studied in the context of carbon nanotube [21|l45] . It was shown that 
a cylindrical carbon nanotube is characterized by two integer charges, since 
a two dimensional solid has two independent density modulations. The intu- 
itive, geometric approach [21"|l45] that leads to this result is to cut and roll a 
sheet of graphite into a tube with no defects. This method provides a general 
and intuitive way of constructing closed-but-not-exact differential form on 
curved surface. Essentially the same method is also used in the construction 
of nontrivial fiber bundle structures [13j|35j. 

The remaining of this paper is organized as follows. In Scc.[2]wc introduce 
differential forms as a natural description of one dimensional translational 
order on curved surface. We also discuss cohomology theory and apply it to 
smectic order on a torus. In Sec. [3] wc give a more in-depth discussion of 
various topological aspects for smectic systems on curved substrates, as well 
as explore the relation between elastic deformations and exact differential 
forms. In Sec. [H we discuss the energetic aspect of smectic orders in curved 
space, and calculate the ground state phase diagram of defects-free smectic 
order on a thin torus. In Sec. wc study spherical smectic and discuss two 
types of smectic states, i.e. hedgehog/latitudinal states and spiral states, 
on a sphere. In Sec. [6j we use the formalism of holomorphic forms, which 
we develop in Appendix [A] to address the third type of states for spherical 
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smectics, i.e. quasi-baseball states. Finally in Appendices El |B] [C] [D] we 

present some formalisms and calculation details that we used in the text. 



2 Analytic Description of Smectic Order on Curved Surface 

2.1 Differential Forms and Density Wave Modulation 

We shall only consider smectic ordering on two dimensional surfaces in this 
work. The formalism, however, can be straightforwardly generalized to three 
dimensions. A two dimensional curved surface, i.e. a Riemannian manifold, 
embedded in the three dimensional Euclidean space can be parameterized 
locally, using arbitrary curvilinear coordinates x = {x 1 ,^ 2 }, as a vector 
function r(x 1 ,x 2 ) = r(x) which gives the position of the point (.t 1 ,^ 2 ) in 
the embedding space. A smectic order on the surface can be described by 
one dimensional density modulation Eq. JT]). We have already seen in Sec.Q] 
that the phase field 0(x) is generically not globally defined as a real- valued 
function. Nevertheless, if we choose a reference point x with phase 0(x), then 
the phase at a neighboring point x + Sx can be written as 

6>(x + 6x) ps G(x) + ip a (x)6x a , (3) 

as long as Sx is small comparing to all the other characteristic length scales. 
Here {ip a ,a = 1,2} is the local wave vector for the density wave around x, 
expressed in the generally non-orthogonal coordinate system 5x a . 
In modern differential geometrjH, the object 

ip(x) = ^ a dx a (4) 

is called a differential form of rank one, or 1-form in short. The set {dx 1 , dx 2 } 
constitutes a bases of the vector space of all 1-forms. While the classical 
meanings of dx a are infinitesimal displacements of coordinates, in modern 
differential geometry, they are defined as the dual vectors of the tangent 
vectors at point x, much like wave vectors are the duals to real space vectors. 
Let {e Q = dr/dx a } be the natural basis vectors of the tangent plane at 
point x, an arbitrary tangent vector V can be expanded in terms of these as 
V = V^ep. The differential forms dx a act on a tangent vector V as linear 
operators and yield numbers: 

(dx a ,V p e p )=V a , (5) 

i.e., they simply extract the a-th components of the vector V. A general 
differential form Eq. ([4J then acts on the vector V as 

(^,V)=^ a V a . (6) 

For a given 1-form ip at x, let us now consider all tangent vectors V 
at x that satisfy (ip,"V) = 0. It is clear that all these vectors form a one 



3 For a physici sts-friendly introduction of modern differential geometry, see ref- 
erence I35lll3ll49l . 
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Fig. 3 (Color online) The state of black layers and that of pink layers differ only 
by a shift of constant phase, and therefore are described by the same differential 
form. 

dimensional subspace of the two dimensional tangent space. We can think 
of this subspace as the infinitesimal version of a smectic layer that passes 
through x. Let us now consider a 1-form field defined on the manifold, i.e. 
there is a 1-form ip( x ) defined at every point x. We try to find a curve 
r(s) = r(x 1 (s), x 2 (s)) such that its tangent vector ^| = ^-e Q everywhere 
satisfies 

<^J>=^^ = 0. (7) 

Such a curve can always be found in the neighborhood of x, by solving the 
above first order differential equation (which is usually nonlinear in x a ). 
Therefore a 1-form ip(x) locally defines a curve passing through an arbitrary 
given point x via Eq. ([7]). 

For an arbitrary smectic state on a curved substrate with no disclination, 
we can always choose a direction along which the phase field increases. 
The differential of the phase field ip = dO is therefore always well defined, 
even though the phase field itself may not be. On the other hand, for a 
given 1-form field ip, the phase at a point y is related to that at point x via 

9{y) = 9{x) + r V = 0{x) + f V ip a dx a . (8) 

J x J X 

For the moment we shall assume that this integral is independent of the path 
connecting x with y. We shall discuss what happens when this condition is 
violated below. Hence we can pinpoint the exact location of all smectic layers, 
as soon as we know the phase at one single point x and the 1-form field ip. In 
other words, the phase field is determined by the differential form ip only 
up to an unknown constant 0{x). Tuning of this constant is equivalent to 
translating all smectic layers in a common direction, as illustrated in Fig. [3l 
which can not be captured by the differential form ip. 

1-forms therefore provide a natural local description of smectic ordering in 
arbitrary curved manifold. This fact is actually well known in mathematics — 
mathematicians often visualize a 1-form field ip in three dimensional space 
as a set of onion skins (analogue of smectic layers!), and a tangent vector V 
as a needle with given length. The pair (ip,~V) simply measures the number 
of layers pierced through by the needle. 
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The metric tensor 

g = 9a/3 dx a ® dar 
is defined by the inner product between basis vectors: 

g a p = d a r ■ dpr = e a ■ ep. (9) 

As usual, we can use the matrix g a p, as well as its inverse g a @ to raise and 
lower tensor indices. A contravariant vector, i.e. tangent vector, tp associated 
with the differential form ip can be defined as 

^ = g a ^ a ep=^e . (10) 

We shall however use the same notation ip both for the 1-form and for the 
corresponding tangent vector ip below. This should not cause any confusion 
for a Riemannian manifold. The norm, or magnitude, of the form ip, is defined 

by 

The geometric meaning of \i\)\j1'K is the number of smectic layers in unit dis- 
tance along the layer normal. In other words, 2ir/\ip\ is the normal distance 
between two neighboring smectic layers: 

<-w 

In the purely geometric description of smectic orders, we do not dis- 
tinguish the smectic layer normal and the nematic director field. Both are 
characterized by the same unit 1-form: 

h = h a dx a = 4r^ = \i>\~' l ip a dx a . (13) 

m 

We shall also use the same notations for the corresponding unit tangent 
vector field 

h = h a e n = ^rrtp = -i- ij; a e a . (14) 



Clearly 



\h\ 2 = g a pn a h = g af> n a n p = 1, (15) 



where n a = g a °hf3 
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A B 

Fig. 4 A): The integral of ip along a close loop 7 vanishes if there is no defects 
inside 7. B): If there is a dislocation inside 7, the integral is 2n times the dislocation 
charge. 



2.2 Dislocations 

The exterior derivative, or antisymmetric derivative, of the 1-form ip is a 
differential form of rank two given by: 

p = dip = d a ipf} dx a A dx 13 
= (ditp 2 -d 2 ipi)dx 1 Ada; 2 , (16) 

where 

dx 1 A dx 2 = dx 1 ® dx 2 - dx 2 ® dx 1 = -dx 2 A dx 1 . 

Consider a small loop 7 that can be continuously shrunk to a point 0, 
as shown in Fig. 3J According to the definition of the form ip, the contour 
integral of ip along 7 is the total change of the phase field, AO, along 7. 
Using Stoke's theorem we can transform this integral into a two dimensional 
integral over the enclosed region D: 

AO = <f V = / dtp. (17) 

J-f J D 

If there is no defects inside 7, as shown in Fig. the overall change of the 
phase field AO must be zero: 

dip = 0. (18) 

By contrast, if there is a dislocation inside the loop 7, as shown in Fig. [4)3, the 
phase change AO along the loop must be the dislocation charge n multiplied 
by 2tt: 

AO = I ip = I dtp = 27rn. (19) 

J j J D 

Eq. (|19p constitute the quantization condition that is to be satisfied by the 
differential form ip, if it describes a smectic order. The exterior differential 
p = dtp must be identified with the dislocation distribution. The loop 7 
and the integer n are the analogues of Burgers' circuit and Burgers' vector 
in the case of crystalline dislocations. If there are k dislocations at places 



4 Note that this would not be possible if there is hole inside 7. 
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{xi, . . . , Xfe} with charges {ni, . . . , nk}, the total dislocation density is given 

by 

k 

p = dip = 2ir n a <5 2 (x — x a ) dx 1 A <ix 2 , (20) 

a=l 

/"# = 27r^n a . (21) 

a=l 

In modern differential geometry, a differential form -0 satisfying = 
is said to be closed. Hence defects-free smectic states are described by closed 
1-forms. A 1-form ip is said to be exact if it can be written as the differential 
of a function 0: ip = dO. An exact form is clearly closed, since 

dip = d(dO) = {d a df}0 - df}d a O)dx a A dx (i = 0. (22) 

This result is called the Poincare lemma. In the context of smectic ordering, 
it means that a smectic state where the phase field <9(x) is globally well 
defined does not contain any the dislocation, a rather obvious result. 



2.3 Global Dislocations and Cohomology Theory 

In a flat space, if there is no dislocation everywhere, dip = 0, one can find a 
globally defined phase field <9(x) such that ip = dO. Consequently the integral 
Eq. (|17p along arbitrary loop 7 always vanishes. This is the inverse of the 
Poincare lemma, which holds in flat space, but may fail in space with different 
topology. One simple counter example on a cylinder is already shown Fig.[T]of 
Sec. Q] Here we give a more detailed discussion in terms of differential forms. 
Using the cylindrical coordinate system (z, ip), let us consider a smectic state 

ip = a dz + k dtp. (23) 

One easily check that dip = 0, hence this state is dislocation- free. The metric 
tensor Eq. ^ is given by 

g vv = R 2 , g zz = 1 g V z = g zv > = 0. (24) 
The layer spacing can be calculated using Eq. (fTTj) and Eq. (fT2|) : 

• 2 * (25, 



Va 2 + R~ 2 k 2 

which is independent of z and (p. Hence Eq. (f2"3"| indeed describes a defects- 
free smectic state with constant layer spacing. 

The integral of ip over the loop 7 y wrapping around the cylinder once is 
given by 

ip= kdp = 2nk,. (26) 
i v Jo 



12 



Hence k must be an integer; it is the global dislocation charge associated with 
the non-retractable loop j v . As discussed in Sec. [IJ k is also the minimal 
number of smectic layers intersecting the non-retractable loop 7^.0 

The integral Eq. (|26|) is independent of the shape and position of the loop 
"f v , as long as it winds around the cylinder once. To see this, we note that 
the integer k can only jump discretely. However, in the absence of disloca- 
tion, such a jump can not happen as one continuously deforms the contour. 
Likewise, one may continuously deform the smectic layers and the integer k 
does not change. The global dislocation charge k therefore is invariant under 
arbitrary continuous deformation of the smectic state, i.e. it is a topological 
invariant. It classifies all topologically distinct defects- free smectic states on 
a cylinder. 

Global dislocations are consequences of the nontrivial topology of the 
substrate, and are naturally addressed by deRham's cohomology theory [35l 
H31I49] . In a compact manifold, a non-retractable loop is called a 1-cycle, while 
a closed but not exact 1-form is called a 1-cocycle. De Rham's theorem of 
cohomology says that cycles and cocycles are dual to each other, and that 
there are equal number b\ of independent cycles and independent cocycles. 
This integer b\ is called the first Bctti number and is a topological property 
of the manifold. Let {/3i,/?2, • • • , A>i} be the b\ linearly independent cocycles, 
we can always find the b\ cycles {71,72, • ■ ■ , 7&i} that form the dual basis, 
such that 



An arbitrary closed form ip can be expressed as a linear superposition of (3i 
and an exact form d<P: 



As will be discussed in detail in Sec. the exact form d<P describes elas- 
tic deformation, and therefore does not carry any topological content. The 
closed form ip Eq. ([28]) is exact if and only if all the coefficients Ni arc zero. 
Furthermore, if ip describe a smectic order, it must satisfy the quantization 
condition Eq. (|19p . hence all the coefficients Ni in Eq. \28\) must be integers. 
They are precisely the global dislocation charges that classify all topologi- 
cally distinct smectic states. The first Betti number of cylinder is one; hence 
only one nontrivial cycle, i.e., the loop j v winding around the cylinder once, 
and only one nontrivial cocycle, i.e., ip = dip. The 1-form Eq. (|23[) is the 
linear superposition of this cocycle dip and an exact form adz. 

5 Even though the form ip in Eq. ([23} looks like an exact form, it is really not. 
This is because the angular coordinate ip is not a single-valued function on the 
cylinder. Hence there is not a single valued phase field such that ip = d&. The 
integral of an exact form dO over an arbitrary loop 7 vanishes identically, regardless 
of whether 7 is retractable or not. 

6 The de Rham's theorem works for differential forms with arbitrary ranks. What 
we invoke here is the special case of 1-forms. 




(27) 



ijj = Ni fii + . . . + N bl p bl + d$. 



(28) 
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Fig. 5 A): Two non-retractabe loops on a torus. B) Coordinate system for torus. 

-7T < 9 < TV, < (f> < 2-K. 

2.4 Smectic Order on Torus 

Let us apply the cohomology theory to defects-free smectic packing on a 
torus. The first Betti number of a torus is b\ = 2 35, 13,, hence there are two 
independent 1-cycles, 70 and 7^, wrapping around the torus in two different 
directions, as shown in Fig.JSK- Correspondingly there are two independent 
1-cocycles, naturally chosen as dO and d(f>, where 6 and <p are the angular 
coordinates wrapping around the loops 70 and 7^, see Fig. 03. Wc easily see 
that {70,70} is indeed the dual basis of {d8,dtp}: 



dO and dip are not exact, since and <p are not single valued function. 

An arbitrary defects-free smectic state ip on a torus is therefore charac- 
terized two global dislocation charges Ng and N^: 



Two integers (Ng,N,p) classify all topologically distinct, defects- free smectic 
states on a torus, and are the minimal numbers of smectic layers intersecting 
the non-retractable loops 70 and 70 respectively. Two states with different 
charges Ng and JV^ are topologically different and therefore are mutually 
inaccessible by elastic deformations. Three example states are shown in Fig. [6] 
Since ip and —ip describe the same smectic state, according to Eq. (|29|) . 
(Ng 7 N^) and (— ATg, — AT^) correspond to the same topological class of states. 
We shall therefore choose the convention > to avoid double counting. 
Hence a defects-free state Eq. (|29]l can be represented as a point in the lattice 
in the upper half (A^, A^) plane, as shown in Fig. [JJ On the other hand, under 




Ng dd + N lf ,d<j) + d$(9, (j)) 

(Ng + dg$) dO + (N4 + d^) dxj>, 



(29) 



where 
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ABC 

Fig. 6 Smectic packing state on torus with (N g = 0, = 28), (Ng = 12, N$ = 0), 
and (TVg = 11, iV^ = 5) respectively. The last one is left handed. 



right handed 



left handed 




Fig. 7 Smectic states on a torus with distinct global dislocation charges are shown 
as grid points in the upper-half (N$,N$) plane. T he c urve (ellipse) is the locus of 
zero nonlinear strain w = 0, as determined by Eq. (|63[1 , States in the shaded region 
have low strain energy and are approximately degenerate in the limit of large system 
size. States well inside the ellipse are expected to exhibit buckling instability. 

spatial inversion of the three dimensional embedding space, the coordinates 
(0, (j>) transform as 

6->-9, + (30) 



The 1-form Eq. 



then transforms as 



I i>N e ,Ns = -Ng d6 + N <f) d(j> + d<P(-0, 4> + TT). 



(31) 



Since translation of <p by 7r amounts to rotation of the torus by 7r around 
the z axis, we find that a state (Ng,^) is transformed into (—N$,N$) by 
the space inversion. Consequently states with Ng ^ and N^ ^ are not 
invariant under spatial inversion, i.e. they are chiral. We shall call states with 
Ng > left handed, while states with Ng < right handed. Clearly a left 
handed state is transformed into a right hand state under spatial inversion, 
and vice versa. 

To understand the geometric properties of the states Eq. ([29| let us cal- 
culate the layer spacing for the special case with <P = 0. With the coordinate 
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system shown in Fig. [5)3, the torus is parameterized by two angles (9, <f>) as 

((R<j, + Re cos 9) cos cp \ 
(Rf + R e cos 9) sin 4> . (32) 
R g sm9 ) 

The corresponding metric tensor can be easily calculated using Eq. (|9|): 

2 



gee = Re, 



(R^ + Re cos0) 2 , 

94>e = 0. (33) 



We shall also define the aspect ratio 



Re_ 

R<b 



(34) 



Note that < r < 1 if the cylinder does not self- intersect. We shall be 
especially interested in the thin torus limit r <C 1, where r provides a natural 
small parameter that greatly simplifies the analysis. 
Defining two parameters ne and as follows 

ne ^2^Re Ne > "* S 2^* (35) 
the layer spacing of the state Eq. (|29|) can then be calculated using Eq. (fl~2f : 

! + - n * 1 (36) 

(l + rcos6>) 2 / v ; 




2r cos 9 
1 + n 2 ,/^ 



0(r 2 



where we have expanded in terms of the aspect ratio r up to the linear 
order. To the zero-th order therefore the layer spacing is independent of 
coordinates (0,(j>). Hence in the thin torus limit r — > 0, the layer spacing 
is strictly constant for arbitrary charges (Ng,N^). For nonzero aspect ratio 
r, the variation of layer spacing is relatively small if either of the following 
conditions is satisfied: 

r < 1 or n 2 e > n\ . (37) 

These are the cases where we expect that the defects-free smectic states 
Eq. (|2"5)l have low free-energy. By contrast, on a fat torus and if rig <C rA, 
the state Eq. ([2"5|) has significant variation in layer spacing, and therefore 
costs large amount of strain energy. Dislocations are expected to proliferate 
to reduces the elastic strain energy. More detailed study of the energetic 
issues for torroidal smectic will be presented in Sec. |U 

Let us determine the loci of smectic layers for a state Eq. (|29| with 
<P = 0. Without loss of generality, let us assume Ng ^ 0. Let us choose the 
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coordinate system (9,(f>) such that one layer, parameterized by (0(t), (f>(t)), 
passes through point 9 = <fi = 0. This layer then must satisfy the equation: 

whose solution is simply given by 

N e 9 + = 0. (39) 

If one follows this layer and wraps around the torus in 8 direction once, 
one obtains AO = Ng x 2ir, and therefore ends up with Ng-th layer. The 
Ng — 1 layers between 0-th layer and Ng-th. layerl are described by the set 
of equations 

Ng9 + N^ct) = 2kiT, k = l,2,...,N e -l. (40) 

This is in fact how the states in Fig. [5] are generated. 

The formalism of differential forms also allows us to deduce interesting 
properties of smectic states with dislocations on a torus. For example, let i/j 
describe a smectic state on a torus. The total dislocation charge is then given 
by the integral of the dislocation density dtp over the whole manifold: 

iV tot = / # = I yj = 0. (41) 

JT 2 JdT 2 

This is because a compact manifold such as T 2 by definition docs not have 
a boundary. Hence the right hand side of Eq. (|4Tj) must vanish. Therefore 
the total dislocation charge of smectic on a compact surface must be zero. 
Obviously, the same result holds for crystalline order on arbitrary compact 
manifold as well. 



3 What do we mean by topological properties? 

The precise definition of topology as a mathematical concept is quite in- 
volved. Loosely speaking, however, a topological property is invariant un- 
der continuous deformations. Mathematicians are primarily interested in the 
topological property of manifolds themselves, hence they consider continuous 
deformation of the manifolds. Physicists, however, are usually interested in 
the topological properties of field configurations defined on manifolds, the 
relevant continuous deformations are therefore those of the fields on a fixed 
manifold. This distinction is seldom made clear in the literature. Of course, 
topology of fields is intrinsically related to that of the manifold where it is 
defined; by studying one we can learn about the other. This relation is indeed 
the origin of the global dislocation charges discussed in this work. 

We are therefore interested in the properties of low energy smectic states 
on curved surfaces that are invariant under arbitrary continuous deforma- 
tions. Using the purely geometric description, a smectic state described as 
a stack of infinitesimally thin layers. By a continuous deformation we mean 
a continuous distortion of the layers such that none of them gets broken; 
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ABC 




D E F 

Fig. 8 Topology of smectic in flat space. A): ground state; B) defects-free deformed 
state that is topologically identical to A); C) a state with isolated dislocation with 
charge one; D) two dislocation with internal distance two; E) two dislocations with 
internal distance one; F) two dislocation with internal distance two. The deforma- 
tion brings D) to E) is plastic; the one brings D) to F) is plastic. States D) and F) 
are topologically identical. 

such a deformation shall also be called an elastic deformation for obvious 
reasons. We shall define two smectic states to be topologically identical, or 
equivalent, to each other, if there exists an elastic deformation that brings 
one state to the other. By definition elastic deformations do not change any 
topological property. By contrast, any deformation that breaks some layers 
shall be called plastic, and necessarily change some topological properties. 
The concepts of elastic as well as plastic deformation obviously apply to all 
systems with translational orders, and reduce to the conventional definitions 
in the case of crystals. 

A number of results follow directly from the above definitions. Layer spac- 
ing is not a topological quantity; neither is a kink, i.e. slope discontinuity, of 
layers. In a flat space with no boundary, all defects-free states (Fig. [8)3) are 
topologically identical to the perfectly ordered ground state (Fig. [SJA) where 
all layers are flat and equal-spaced. Any state with a single dislocation with 
charge k (Fig. [Sp) is clearly topologically distinct from the ground state. 
Two such states with distinct charges are also topologically distinct from 
each other. The dislocation charge is therefore a topological property. For 
states with two dislocations (Fig. [5p, E, F), the number of layers between 
them is a topological property. This number is the "internal" distance Q, i.e. 
the distance measured by the local smectic layer spacing, between them. To 
change this internal distance (Fig. [Hp to E), it is necessary to break and 
reconnect smectic layers. In contrast, the physical distance between dislo- 
cations can be changed by changing layer spacing and without change of 
the internal distance (from Fig.[8]D to F); such a deformation is elastic. For 



7 We adopt this terminology from Kroner's study [26] of crystalline defects using 
differential geometry. 
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states with many dislocations, the internal distance between two dislocations 
depends on the path connecting them. 

A plastic deformation that breaks infinite number of layers shall be called 
a non-local plastic deformation. For example, one only needs to cut three 
layers in order to go from the defects free state Fig. [SJ'V to the state Fig. [Sp 
with a pair of opposite dislocation pair. Such a process is local, and can 
spontaneously happen in thermal equilibrium. On the other hand, it takes a 
nonlocal plastic deformation to go from Fig.[8j\. to Fig.[8fU: one needs to cut 
infinite number of layers in order to introduce an isolated dislocation into a 
defects free state. 

The homotopy theory of defects [33" l l34 p 23 p 24| . developed in 1970's, deals 
with deals with defects in various ordered media in Euclidean spaces. In this 
theory, a topological property is defined to be invariant under arbitrary local 
deformations. For systems with translational orders, such as smectic or crys- 
tals, these local deformations include both elastic and plastic deformations. 
For example, the states in Fig. [8JA.&D are topologically equivalent according 
to the homotopy theory. Topological properties according to our definition 
are therefore not necessarily topological according to the homotopy theory. 

As stressed by Mermin [33J long ago, the homotopy theory of topologi- 
cal defects, though rigorous and elegent, yields lots of strange results when 
applied to translational orders. For example, by transporting a dislocation 
around a disclination, the Burgers vector of the former can be changed. Con- 
sequently, dislocation almost completely loses its topological identity in the 
presence of disclinations. The Homotopy theory therefore is not able to distin- 
guish states with same disclination charge, but otherwise are widely different. 
While these results are not A priori wrong, they are very counter-intuitive 
and often misleading. It is therefore fair to say that the homotopy theory 
is not refined enough to capture the differences and connections between 
oricntational defects and the translational ones in a translationally ordered 
system. This point may become clearer when we discuss smectic order on 
sphere in Scc.[5]t[ni 

Differential forms, augmented by de Rham's cohomology theory, as dis- 
cussed in this work, provides a natural frameworks for study of translational 
orders both in flat space and in curved space. The topological significance of 
defects pairs, as shown in Fig. [8)3 , E &F, is captured by differential forms 
that are not closed; elastic deformations are naturally described by exact dif- 
ferential forms. Being a geometric description, differential forms also provide 
a natural description for the energetics of translational ordered systems. The 
issue of disclinations is more subtle, and in principle should be described by 
nontrivial fiber bundle structure associated with the differential forms. For 
two dimensional smectics, however, there is an easy alternative in terms of 
multiple valued complex differential forms, as will be discussed in Sec. [6] 
A systematic exploration of disclinations in translational systems in terms 
of fiber bundle theory shall be left to a future work. Finally we note that 
this geometry-based approach is in the same vein as the intuitive Volterra 
treatment of dislocation and disclination in continuum solids. 



For system with orientational orders, the distinction between plastic deforma- 
tion and elastic deformation is irrelevant. 
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Fig. 9 (Color online) Yellow lines are undeformed smectic layers. Red lines are de- 
formed layers. The point (dot) with coordinates x is deformed to the point (square) 
with coordinates x. x and x are therefore the Eulerian as well as Lagrangian coor- 
dinates for the deformed point (square). 

3.1 Elastic Deformations as Exact Differential Forms 

It is a well known result in differential geometry that the integral of a closed 
1-form along a loop is a topological quantity, that is, it is invariant under 
arbitrary continuous deformation of the loop and of the differential form 
itself. This is the mathematical origin of the global dislocation charges that 
we discuss in the preceding section. □ In this section, we study in detail how 
elastic deformations change the differential form describing the smectic order. 

Let us consider a defects-free smectic state on a torus, described by a dif- 
ferential form Eq. (|29)) . It is characterized by two global dislocation charges 
(Ne,N$). Let us elastically deform the substrate, together with the smectic 
pattern, such that a point with (9, </>) is brought to (9, (f>) by the deformation, 
as illustrated schematically in Fig. [5J Using the convention of elasticity the- 
ory, we shall call (9, (f>) the Eulerian coordinates and (9, </>) the Lagrangian 
coordinates. The elastic deformation then define (9, <j>) as functions of (9, <f>) 
and vice versa. The elastically deformed smectic state can then be obtained 
by simply replacing in Eq. ([25]) the coordinates (9, 4>) by the Lagrangian co- 
ordinates (0,0), which shall in turn be treated as function of the Eulerian 
coordinates {9,4>): 

V>def = N e d9{9, 0) + dj>(9, <j>) + d${9, $) (42) 

Let us now look at the difference Aip between ip^ei and i/j: 

Alp = V>dcf - "0 

= N g (d9 - dff) + N^d} - d<jy) 

+ d ($(6, 4>) - $(9, <p)) . (43) 

9 Similar phenomena also arise in lattice gauge theory [25] with periodic boundary 
conditions. The relevant space is then the d-dimensional torus, with d independent 
non-retractable loops. The integral of the gauge potential (a 1-form) along these 
loops then give d topological charges that are invariant under arbitrary gauge trans- 
formations. Hence the configuration space consists of many components that are 
not mutual accessible by any local operation. 
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Because the Eulerian coordinates and the Langrangian coordinates are re- 
lated by elastic deformation, the angles 9 (</>) wraps around the torus if and 
only if 9 (<j>) does the same. This implies that the integrals of Aip both along 
the non-retractable loop jg and along 7^ vanish. Hence, according to the de 
Rham's theorem, the integral of Aip along arbitrary loop vanish and Aip is 
an exact differential form. Therefore the effect of an elastic deformation is 
to add an exact 1-form Aip to the original differential form ip, and does not 
change either global dislocation charge. 

An infinitesimal elastic deformation can be generated by a vector field 
X = X a e a . Here for convenience we use the shorthand x 1 = 9, x 2 = (p. 
Mathematically this means that a point x (Lagrangian coordinates) is de- 
formed to be the point x (Eulerian coordinates) such that 

x a =x a + eX a , (44) 

where e is the infinitesimal flow parameter, and X a arc considered as func- 
tions of the Lagrangian coordinates x 13 . The deformed 1-form is then given 

by 

</, dof = ip a ( x )dx a = ip a (x - eX)d(x a -eX a ) 

= ijj a (x)dx a - e {ip a d p X a + X a d ip a ) dx 13 , (45) 

where we have only kept terms up to linear order in e. We recognize that the 
coefficient of e is negative of the Lie derivative of the 1-form ip. Hence the 
infinitesimal change of ip induced by the elastic deformation is given by 

Aip — ip dei - ip = -eCy^ip. (46) 

Using the infinitesimal homotopy relation [35 l lT3 l l49] 

£x = dix + i x d, (47) 



where ix is the interior product, as well as the fact that dip = 0, Eq. (|46|) 
reduces to 

Aip = d (-e i x ip) = d (-e ip a X a ) . (48) 

We have therefore explicitly constructed the potential corresponding to the 
exact form Aip that arises from an elastic deformation generated by the 
vector field X. 

Conversely, an arbitrary change of the 1-form ip by an exact form d<P can 
be generated by an elastic deformation. This elastic deformation is however 
not unique. In the case of torus, for example we can write Eq. (|29[) as 

= N e d \6 + -^${0, <f>)] + N^dcP 

= N e d§ + d<f>, (49) 

where 

= 6+ ±$(0,4), ^ = cp. (50) 
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According to our preceding discussion, then, this smetic state tp can be gen- 
erated from the prototypical state 

V'orignal 

NgdO + N^dcj) (51) 

by the elastic deformation 

0,$)—>(e,<t>), (52) 

Again (9, <f>) are the Eulerian coordinates and while (9, <p) are the Lagrangian 
coordinates. This result provides a starting point for systematic studies of 
thermal elastic fluctuations in defects-free, translationally ordered system on 
a torus. 



4 Geometry and Energetics of Smectic on Curved Substrates 

Our problem is to study the topological as well as geometric properties of the 
low energy smectic states on a rigid and compact substrate, This mission can 
be decomposed into two steps. Firstly we use de Rham's theorem to classify 
all low energy smectic states with minimal number of defects. For the case 
of torus, for example, this is achieved by Eq (|29|) . Note that the first step 
analysis is independent of energetics. As the second step, we minimize the 
total free energy, for given values of the global dislocation charges Ni , over the 
exact form d<P. Physically, this means we let the system elastically relaxed to 
the energy minimum in a given topological sector. A similar two-step process 
is also taken in the linear theory of topological defects in flat space, for 
example, in smectic and in crystals [5ll37| , where one first finds a singular part 
of the phonon field that gives the right topological charge, then adds a smooth 
part to this phonon field to obtain the energy-minimizing solution. The same 
strategy is also implicitly employed in the Volterra's intuitive construction 
of defects in solids. 

To carry out the second step, however, we do need to understand the 
energetics of smectic on an arbitrary rigid substrate. The total elastic free 
energy of a smectic consists of two parts. Firstly deviation of the layer spacing 
away from its energetically preferred value do costs elastic free energy. The 
layer compression/dialation can be captured by a dimensionless nonlinear 
strain w defined as 

«>(*)=(l-|) = (l-^M a ), (53) 

where used was Eq. (|12[) . The sign of w is chosen according to the con- 
vention in elasticity theory, i.e. w > (w < 0) corresponds to a dilational 
(compressional) strain, where the layer spacing increases (decreases). The 
corresponding strain free energy is 



Fb = ^ dcjw(x) 2 , 



(54) 
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where dio = yjg dx 1 Adx 2 is the invariant volume form. It can be easily shown 
that this strain is identical to the one used in reference [47] , and also reduces 
to the well known nonlinear strain for smectic order in flat space p~5ll8] . 

Secondly, variation in the layer normal, or bending of the layers, also 
costs elastic free energy. The elastic free energy of vector/nematic/hcxatic 
orders constrained to the tangent space of a flexible membrane has been 
explored extensively by many authors [111130.22.39. 12]. The most general 
form of this elastic free energy contains many independent parameters. Here 
we shall focus on the simpler case of smectic orders on a rigid substrate, 
and leave the more complicated case of flexible substrate for a future work. 
We first emphasize that it is the variation of the layer normal in the three 
dimensional physical space that costs elastic free energy. In general, the elastic 
free energy depends on the director field n as well as its covariant derivatives 
D a n. It can also depend on the metric tensor g a p, the antisymmetric tensor 
e a ' 3 = \/<? _1 7°^, with 7 Q/3 the antisymmetric symbol, as well as the the 
extrinsic curvature tensor K with components K a p = N ■ D Q e^, where 

N=^^r, (55) 
|ei x e 2 \ 

is the unit normal vector of the surface. In the present case of a rigid sub- 
strate, however, these latter objects are just fixed parameters, not dynamic 
quantities. All tensor indices must be properly contracted over to ensure the 
reparameterization invariance. Furthermore, the tt rotation symmetry per- 
taining to the nematic order requires that the elastic free energy must be 
invariant under the transformation h — > — n. This is in strong contrast with 
the tilt order in double layers extensively studied long time ago [39]. Fur- 
thermore, we shall consider only surfaces whose two sides are equivalent, i.e., 
those invariant under N — > —N, so that terms like h ■ K • n are not allowed 
in the free energy. 

Incorporating all these constraints we find that there are six terms in 
total which arc quadratic in h and are marginally relevant under naive power 
counting^- First of all there are two terms involving the covariant derivatives 
of the director field: 



Fi J = dw 



(56) 



where K\ and K% are splay and (inplane) bending Frank constants. Note 
that 

(D x h) = e al3 D a hp (57) 

is a pseudo-scalar in two dimensions. These two terms can be obtained by 
naively promoting the 2D flat space Frank free energy for nematics to curved 
space straightforwardly. In the absence of dislocation, however, the inclusion 
of bending term with coefficient K3 is not essential, since as is well known 

10 We note that these terms can also be expressed in terms the unit tangent vector 
t of the smectic layers, where t • h = 0. 
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(also see discussion in Sec. [SJ, bending of nematic director necessarily leads 
to change of smectic layer spacing, which is already penalized by the strain 
energy Eq. ([54]) [46] . 

Secondly there are three independent terms quadratic in the extrinsic 
curvature tensor K af3 [22l[TTll39lfT2] : 

y (n ■ K • hf + J 2 {h • K ■ t) 2 + ^-(t • K • tf, (58) 

where t is the unit vector tangent to the smectic layers, i.e. perpendicular 
to n, with i a = e£iip, while raising and lowering of indices using the metric 
tensor is implicit in the above expressions. It is interesting to note that in 
the special case of a sphere, the extrinsic curvature tensor is isotropic, hence 
all these three terms reduce to trivial constants and therefore can be safely 
ignored. 

Finally for chiral systems, a term linear in e a p is also allowed: 

ch ■ e ■ K • n — ct ■ K ■ n. (59) 

We shall, however, not consider the effect of chiralty in this work. 

All the above terms were explicitly listed in reference [3S] for the case 
of symmetric bi-layer with tilt order. Also discussed there are some other 
terms that are allowed for tilt order, but not for nematic order. For hcxatic 
membranes, none of the terms in Eq. (|58|) or in Eq. (|59p is allowed. Further- 
more, the sixfold rotation symmetry of hexatic order requires K\ and K$ in 
Eq. ([5"6"]l to be the same. 

Geometrically, the terms in Eq. ([58]) express nothing but the energy cost 
for the director bending due to the extrinsic curvature of the substrate. There 
is no symmetry principle that allows us to determine the relative importance 
of coefficients of the above terms E3- Nevertheless, thermodynamic stability 
does require both K\ and K 3 in Eq. ([56]) to be positive. Otherwise, the 
system would develop in-plane nematic director modulation spontaneously. 
The coefficients Ji in Eq. ([58]) are more subtle: if J\ = J2 = Jz ~ J, we easily 
see that Eq. ([58]) reduces to \ JTrK-K, which is independent of the dynamic 
variables n and t. For a fluctuating liquid crystalline film, this term simply 
rcnormalizcs the clastic modulus of extrinsic curvature. By contrast, for a 
rigid substrate, which is our primary interest in this work, this term is just an 
irrelevant constant that can be safely ignored. For the latter case, therefore, 
only the differences between Ji 's have physical meanings and the signs of the 
coefficients Ji can not be determined by any physical consideration. 

Below we shall consider a particularly simple limit where J\ — 2 J2 — 
K\ = K3 = K and J3 = 0. Simple calculation then shows that Eq. ([58]) 
reduces to K/2 h ■ K 2 • h, and the sum of Eq. ([56]) and Eq. ([58]) reduces to: 

Fn = y J du g a ^d a n ■ d n, (60) 

1 The elastic model considered in [47] corresponds to K\ = J3 = K and all other 
moduli vanish. A more general form of elastic model is considered in [12] . 
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where d a are the ordinary derivatives, while n is constrained to be tangential 
to the substrate. We note that Eq. ([60|) is the elastic free energy for a nonlinear 
Sigma model in a two dimensional curved space, with the constrain that the 
director n must be tangential to the substrate. For K > 0, Eq. (|60[) penalize 
the bending of the director n in the three dimensional physical space. As a 
consequence, the nematic director would prefer to align with the flatter axis 
of principle radius, with larger radius of curvature. This is indeed seen in 
recent experiments on smectic vesicle [55]. By contrast, in the model studied 
by Santangelo et. al. [47] the smectic layers prefer to align along the flatter 
direction. These two models clearly describe different systems. 

The total elastic free energy for smectic on curved substrates is then 
given by the sum of Eq. ([M]) and Eq. ([6H| . characterized by two independent 
parameters B and K . 

4.1 Energetics of Toroidal Smectics 

Let us calculate the elastic free energy of all defects-free smectic states 
(Ng,N$) on a torus. Strictly speaking, we have to calculate the total free 
energy for the generic state Eq. (|29|) with arbitrary <P and minimize the total 
free energy over the elastic deformation d<P(6, </>) for given charges (Ng, N^,). 
To simplify analysis, however, we shall instead only consider states with <P 
set to zero. As we have seen in Eq. (|36j) . in the regime of thin torus r< 1, 
these states have almost constant layer spacing. Our results are not expected 
to be changed qualitatively by including <P into the calculation. 

The nonlinear strain Eq. ([53]) can be easily calculated using Eq. ([36]) : 

w = l-dl(n 2 e + - " ) (61) 

u I (1 + rcos6>) 2 J v ' 

Let us first look at the infinitely thin torus limit r — > 0, where the strain 
becomes independent of 9 and <f>. The elastic strain energy then reduces to 

F B = J cLo^w 2 -iBA(dgn2 + dgn$-l) a , (62) 

where A 4tt 2 RqR^, is the total area of the torus. 

For arbitrary given values of Re and if Ng,N^, were continuous vari- 
ables instead of integers, we would be able to find a one-parameter family of 
states with vanishing strain: 

ivq = a%n 2 B + din 2 - 1 

dlN 2 d 2 N 2 
{2nR e ) 2 (2W? ) 2 
= 0. (63) 

These "states" trace out an ellipse in the upper- half (Ng,N$) plane F 2 !. as 
illustrated in Fig. [7] In reality, of course, (Ng,N^,) can only take integer 



Remember that we choose the convention iV^ > 0. 
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values. For generic Rg and R^, there is no integer solution to Eq. ([63]) for 
(Ng,N^), hence smectic state with strictly zero strain gencrically does not 
exist. There are, however, many low energy states where the strain energy 
is negligibly small. Consider the region in the (Ng,N<p) plane between two 
ellipses: 

44 + 44 = i±-^, ( 64 ) 

which is illustrated as the shaded region in Fig. [71 For all states inside this 
region, the absolute value of nonlinear strain is bounded by do/VA, hence 
the total strain energy Eq. (|62p is bounded by a microscopic value: 

F b <\b4, (65) 

which is negligible comparing with all other free energy contributions, espe- 
cially in the large system size regime. The total number of these low strain 
energy states, i.e. the total number of grid points inside the shaded region 
in Fig. [71 is roughly the total area of this low-strain-energy region in the 
upper-half {Ng^N^) plane, which can be easily calculated: 

M{F b <\b4) = 2^. (66) 

Hence the total number of low strain energy states scales as the square root 
of the torus area. 

Hereafter we shall restrict ourself in the shaded region with low strain 
energy. Let us now consider the effect of a finite but small aspect ratio r. 
The strain energy up to the order of r 2 can be calculated using Eq. ([6"Tj) : 

F B = l -BA (wl + r 2 d 2 nl(2d 2 nl - w )) , (67) 

where the reference strain wq is defined in Eq. f|63[) . 

We have already shown that inside the shaded region in Fig. ([7]), |wo| < 
d /y/A <C 1, and therefore can be ignored in Eq. ([6"7j) . Ej Using Eq. ([6"3"]) . 
then, we can parameterize ng and n^, in the following way: 

do ng « sin a 
d n^ « cos a 

It is not difficult to see that the geometric meaning of a is the angle be- 
tween the nematic director n and the local (f> axis in the thin torus limit, or 
equivalently, the angle between the smectic layers and the local 6 axis. 

13 Need to show that the buckling instability does not set in to the shaded region 
as long as the do < 2£, where £ = y/K/B is the nematic penetration length. This 
condition should always be satisfied. [10] , 

1 Its inclusion, however, does not qualitatively change our results below. 



ng 
nl + n 2 . 



nl + n 2 
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The elastic strain energy Eq. (|67|) can then be simplified to 

-Fb = BA r 2 (cos 4 a + 0(r)) . (68) 

As expected, the strain energy density scales as r 2 and becomes small for 
a thin torus. It is in this limit that states with different global dislocation 
charges, as shown in the shaded region in Fig. [71 are approximately degen- 
erate. Finally, the strain energy vanishes at a = ±7r/2, i.e. all smcctic layers 
are along <j> axis, and the smectic layers spacing are equal everywhere. 

The nematic distortion energy can be similarly calculated, using Eq. (|60[) 
and Eq. (TH)) . Again expanding to the leading order of the aspect ratio r, we 
find 

Fa = (sin 2 a + O(r)) , (69) 

which clearly prefers a = 0, 7r, for which the nematic director is parallel to 
the </> direction with smaller curvature. Note, again, that Fa depends on the 
aspect ratio r but not on the substrate area A, hence becomes less important 
for larger system. 

Let us introduce a dimcnsionless parameter 



where £ = JKjB is usually called the nematic penetration length and is 
of order of layer spacing do in the strongly anchoring limit. The total free 
energy, being the sum of the strain energy Eq. (|68p and the Frank free energy 
Eq. (|69)) . is given by 

f 2 

-F tot = cos 4 a + — sin 2 a, (71) 



r 2 BA lox x 2 

= l + 2(i-l)a 2 + (| * )</ + ..., 

where we have also expanded the total free energy in terms of small angle a. 
Note that the quadratic coefficient changes sign at \ — 1- 

The ground state of smectic order on a thin torus is therefore determined 
by the competition between the strain energy and the Frank free energy. 
Minimization of the total free energy Eq. (fTTj) over the angle a is straight- 
forward. We find two different phases. When % < 1, the Frank free energy 
dominates the strain energy and the total free energy has only one minimum 
at a = 0. As illustrated in Fig. [6^, this corresponds to an achiral phase with 
all smectic layers forming loops in 9 direction. The nematic director is an- 
chored by the extrinsic curvature to the local flatter <f> direction everywhere. 
In this phase, the total free energy is given by BAr 2 , predominantly from 
the variation of layer spacing. On the other hand, when % > 1, the total free 
energy has two degenerate minima at 

a = ± arccosx -1 - (72) 
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Fig. 10 The phase diagram of smectic on a thin torus. The phase boundary is 
given by Rg/R^ = 1. 

As illustrated in Fig. OU, this corresponds to a chiral phase with the ncmatic 
director continuously varying with the parameter \- The spatial inversion 
symmetry is spontaneously broken in this phase, due to the competition 
between the strain energy and the Frank free energy. Since a continuously 
increase from zero as x increases from unity, the phase boundary is of the 
second order. The phase diagram of a smectic on thin torus is summarized 
in Fig. [TO] 

According to Eq. (|7D|) if we fix the aspect ratio r and let Rg approach 
infinity, we always end up with % > 1, hence the system always falls into the 
chiral phase. On the other hand, if we fix Rg and let R^ approach infinite, 
we always end up with \ < 1, and the system falls into the achiral phase. 
This subtlety in taking the thermodynamic limit originates from the fact 
that the Frank free energy does not scale with the system size in the usual 
way. On the other hand, for fixed system sizes Rg <C i?0, if the smectic order 
is sufficiently weak, the nematic penetration length may be large enough to 
make the parameter x < 1, according to Eq. (|T0|) . Finally, it may be a little 
surprising that, according to Eq. (|72|) and Eq. (|70|) . as long as the Frank 
constant K does not vanish, the angle a is always less than tt/2. That is, the 
achiral state with = can never be the ground state. 

We must, however, be careful with the interpretation of the phase dia- 
gram Fig. [101 States with different a are separated by topological barriers. 
Transitions between these states are therefore very slow. The energetic dif- 
ferences between them scale sublinearly with the system size and therefore 
may not play significant role in physical processes. 



5 Smectic Order on Sphere 

5.1 Spherical Nematics and Disclinations 

The first Bctti number b\ of a sphere is zero. Indeed, there is no non- 
retractable loop on a sphere: any loop on a sphere can be continuously shrunk 
into a single point. Hence one would naively expect that the smectic ground 
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Fig. 11 Generation of a one-parameter-family of degenerate ground states for 
spherical nematics in the limit K3/K1 — > 00. One first cut the sphere along a 
great circle that passes both poles, then rotates the two halves by an arbitrary 
angle. The resulting state has four +1/2 disclinations on a great circle and forming 
a rectangle. Note that the director field is everywhere continuous except at four 
disclination cores. 

state on a sphere is topologically unique. This is, however, incorrect, due 
to the presence of orientational defects, i.e. disclinations. By Eq. (fTT)|) a dif- 
ferential form ijj can also be understood as a tangent vector field i/j on the 
sphere. The Gauss-Bonnct-Poincare theorem dictates that the total disclina- 
tion charge of a tangent vector field on a sphere must equal to two. More 
importantly, ip and — tp corresponds to the same state, hence if> really de- 
fines a nematic order on sphere that admits half integer disclinations. In 
the presence of the latter objects, we will not be able to consistently choose 
a direction along which the local phase field increases E3, and therefore ip 
becomes a multiple- valued differential form. 

The ground state of spherical nematics has been subject to intensive study 
recently [53"Il4"5I [T] . and is expected to contain four +1/2 disclinations in the 
ground state. When the bending constant K3 and splay constant K\ equal to 
each other, four disclinations form a regular tetrahedron in the ground state. 
By contrast, in the limit K3/K1 — > or 00, all bending (splay)dcformations 
are forbidden 0- One can then show [48] that all four disclinations form 
a rectangle with arbitrary shape on a great circle. The system therefore 
is characterized by a one-parameter family of degenerate ground states F^L 
Interestingly, all these states can be generated from the special states with 
two integer disclinations at two poles, by a cut-and-rotate surgery, similar 
in spirit to the Volterra construction of dislocation and disclinations in solid 
materials. This surgery is illustrated in Fig. QT] The manifold of degenerate 
ground states is therefore parameterized by the angle of rotation. 

It is well known [15] that the bending constant K$ (for the nematic di- 
rector) is effectively infinity in the smectic phase in 3D flat space, if the layer 
compression modulus is large and if dislocations are energetically expensive. 
Geometrically it is easy to see that bending of director field either leads to 

15 Remember that ip is the differential of the local phase field of the density wave. 

16 We note that on an arbitrary curved substrate, bending-free (splay-free) ne- 
matic texture may not exist everywhere. The system nevertheless can exhibit do- 
main structures. In each domain the nematic director is bending-free (splay- free) , 
while the director is discontinuous on boundaries between neighboring domains. 

17 This degeneracy is expected to disappear as soon as the ratio K3/K1 deviates 
away from two fixed point, or 00. 
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change of layer spacing, or proliferation of dislocations. The same result also 
holds in two dimensional curved space and can be elegantly derived using 
differential forms. Let us write ip = \ip\h, where n = n a dx a is the 1-form 
associated with the nematic director field, while the magnitude of tp, is 
inversely proportional to smectic layer spacing, according to Eq. (fT2|) . Taking 
the exterior differential of ip, as well as using the Lebiniz rule, we find 



Clearly d\ip\ = (d a \ip\)dx a describes change of smectic layer spacing, while 
dh = {d\fi2 — dznxjdx A dx 2 is a 2-form describing bending deformation 
of the director field. Hence if there is no dislocation, dip — 0, and the layer 
spacing is constant, d\ijj\ = 0, the bending deformation of the director field 
is strictly forbidden, dh = 0, which implies that the bending constant is 
effectively infinity. 

Conversely, starting from an arbitrary bending-free ground state for spher- 
ical nematics, it is intuitively plausible that one can glow, layer by layer, a 
smectic pattern with equal layer spacing and with no dislocation. We will, 
however, have to fine tune the layer spacing d around the energetically pre- 
ferred value do so that the smectic pattern can be fit onto the sphere with 
given radius R. For a large sphere R/dn 3> 1, the resulting dimensionlcss 
strain Eq. ([53]) should scale as do/R Ej, hence the total strain energy is 
bounded by 



The key feature of Eq. (|74[) is that it does not scale with the system size, 
and therefore becomes negligible compared with other energy scale as R/dy 
becomes large. Furthermore, all the smectic states thus constructed natu- 
rally have equal Frank free energy. Therefore, similar to toroidal smectics, a 
smectic order on sphere also exhibits many approximately degenerate ground 
states, with the free energy difference bounded by a microscopic energy scale 
Eq. J7^p . Unlike the case of smectic on torus, we are not able to explicitly 
calculate the total free energy for each smectic state on a sphere. We shall 
calculate the total number of low energy states below. In so far as the nematic 
director is concerned, therefore, a spherical smectic behaves the same as a 
spherical nematics in the limit of K3/K1 — > 00, and is expected to exhibit 
four +1/2 disclinations sitting almost on a great circle o This, however, 
does not mean that these two systems are completely equivalent. As we shall 
show below, the presence of smectic layers imposes extra structure on the 
top of the disclination pattern, which can not be captured by the director 
field alone. 

18 The system may relax this strain, in a way that depends on the nematic texture. 

19 We note further that our analysis is of mean field nature, and completely ig- 
nores long wave length fluctuations at scales comparable with the substrate. These 
fluctuations may further break the approximate degeneracy. A calculation of these 
fluctuations is however considerably difficult and shall not be addressed in this 
work. 

20 This result is probably first pointed out by Kleman and Blanc [4]. 



# = d(\ip\ri) = (d\tp\) A n + \ip\dh. 



(73) 




(74) 
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Fig. 12 Three types of spherical smectic states discovered in recent numerical 
simulations. 

5.2 Spherical Smectic and Ambiguity in Location of Disclinations 

Several numerical studies of spherical smectic orders in di-block copolymer 
films confined on sphere have been recently reported [9"I l2"7Il5Tj] . In agreement 
with our discussion above, all the low energy states are found to have ap- 
proximately equal layer-spacing and are free of isolated dislocation. There 
may be, however, at most one layer termination (dislocation of charge one) 
in the disclination cores. If there are two terminated layers near one discli- 
nation core, it saves energy to connect them. This process is of course not 
elastic according to our definition in Sec. [3l Nevertheless it happens with 
finite probability in the presence of thermal fluctuations, since only local re- 
arrangement of the layer structure is needed. We therefore infer that in all 
low energy smectic states on a sphere, there can be at most two open smectic 
layers, with possibly four layer terminations. All other smectic layers must 
form close loops. 

The low energy states of smectic states are empirically classified into 
three categories in reference [9] : The first class is so-called hedgehog states in 
reference [9] EE shown in Fig. [T2R, where there are two +1 disclinations at 
north pole and south pole respectively and all layers are circles of constant 
latitude. We shall however refer to them as latitudinal states. A latitudinal 
state is free of layer termination everywhere. The second class of state, as 
illustrated in Fig. [12)3, is called spiral states, which have the same disclination 
pattern as the latitudinal states, but all smectic layers are spirals around the 
two poles. A spiral state may have one layer termination in the center of 
each disclination. The third class of state is called quasi-baseball state, as 
illustrated in Fig. [T2lC. as well as Fig. [TBI where four +1/2 disclinations sit 
on a great circle and form a rectangle. Again at most one layer termination 
can appear in the core of each disclination. It is clear that each class of states 
contains many different states that are topologically distinct to each other. 
Finally in reference [9] the total free energy of these three classes of states 
are also calculated; it is find that the free energy difference decreases as the 
system size increases, in agreement with our theoretical result discussed in 
the preceding subsection. 

In the remaining part of this paper, we shall explore the topological as 
well geometric properties of these states, as well as the connection between 



The term hedgehog here is likely to be a misnomer. 
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C D 

Fig. 13 Topologically one +1 disclination is not distinguishable with two +1/2 
disclinations. Hence A and B are topologically identical, while C and D are also 
topologically identical. A and C are however topologically distinct to each other: 
There is no continuous deformation that can take state A to C. On a sphere these 
continuous deformations involve significant changes of layer spacing, and therefore 
incur high elastic energy penalty. 




+ 1 

A B 

Fig. 14 Starting from a latitudinal state (A) with N layers and with two +1 discli- 
nations at poles, we split each integer disclination into a pair of +1/2 disclinations, 
and drag them down (in perpendicular directions) to the equator. The resulting 
state (B) is a quasi-baseball state Ni — N2 = N. Only two +1/2 are visible in (B). 

these states. In particular, we shall show that both latitudinal states and 
spiral states can be understood as the limiting case of quasi-baseball states, 
where one of two dislocation charges becomes small or vanishes identically. 

To understand the topological properties of smectic order on sphere, it is 
important to discuss the role played by disclinations. Disclinations in smectic 
are not completely characterized by their disclination charges. For example, 
states Fig. [T3R and Fig. [T3C are topologically distinct with each other, even 
though both have disclination charge +1: they can not be transformed into 
each other by elastic deformation. In fact, they can be obtained by projecting 
stereographically a latitudinal state and a spiral state of spherical smectic 
onto the complex plane. We shall see later that state the disclination C and D 
should be understood as a composite defects with disclination charge +1 and 
dislocation chargs —1. Composite defects with arbitrary dislocation charges 
k can also be constructed in a straightforward way. 
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It is also important to note that specification of layer positions do not 
completely determine the position of disclinations. In all three states shown in 
Fig. 1121 we can only determine the location of disclinations up to the precision 
of layer spacing. This problem of course also exists in flat space. For example, 
for states A and B shown in Fig. [T3J we are only certain that the total 
disclination charge enclosed by the innermost loop is +1, the precise position 
of the disclination is however can not be determined from the layers. More 
interestingly, by continuously deforming smectic layers from state Fig. I131 \ 
to state Fig. I13B . we may say that we break one +1 disclination into a pair 
of +1/2 disclinations and separate them far apart. There is however, no truly 
topological difference between these two states. Similarly the state Fig. [T3T J 
can also be elastically deformed into Fig.fTST). Return to the case of spherical 
smcctics, such splitting and dragging of disclinations allows us to construct 
a quasi-baseball state from a latitudinal state by continuous deformation 
only. This process is shown in Fig. [14] Hence these two states Fig. [T4R 
and B, though seemingly very different, are really topologically identical to 
each other. As we will discuss in detail below, splitting and displacement 
of disclinations provide the basic elastic mechanism that connects together 
many apparently different states of spherical smectics. Nevertheless, we must 
note that the process of splitting and separating defects may change the layer 
spacing significantly, and therefore costs elastic free energy. For example, the 
layer spacing in the final state shown in Fig. [14] is only half of that of the 
initial state. The corresponding free energy cost is therefore proportional to 
the system area, and is enormous for large system, where R/cIq 3> 1. 

Consider a quasi-baseball state on a big sphere R/do 3> 1, such that 
disclinations are well separated. Displacement of disclinations by large dis- 
tance induces significant change of layer spacing which costs high elastic free 
energy, and therefore is strongly suppressed even in the presence of ther- 
mal fluctuations. Each disclination therefore fluctuates slightly around their 
average position. To simplify the problem, we can treat disclinations as im- 
mobile and consider the subspace of all elastic deformations that do not 
change the position of any disclination, and to classify all low energy states 
with fixed positions of disclinations. Two states are then called weakly topo- 
logically equivalent if they can be reached from each other by any elastic 
deformation that leave all disclinations fixed. With this definition, we shall 
find that quasi-baseball states can also be classified by two integer global 
dislocation charges, much similar to the torus case. This similarity between 
toroidal smectic and spherical smectic (with fixed disclinations) turns out to 
have interesting connection with the theory of elliptic functions, as will be 
shown in Sec. [5] The stability of these states however, are not purely of topo- 
logical nature. Rather, they are stabilized by energetic as well as entropic 
barriers. When large scale motion of disclinations is allowed, many of these 
states actually form a continuum, mutually accessible by elastic deformtions. 

By a strong contrast, if two +1/2 disclinations form a close bound pair, it 
does not cost high energy to move one relative to the other. In particular, we 
can twist a disclination pair so as to exchange their positions. This process 
only induces a very small change of the layer spacing (or order of do /N, where 
N is the total number of layers), and therefore only costs small amount 
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of elastic free energy. Swap of disclination pairs therefore are low energy 
excitations in these states and should be able to be observed in simulations 
as well as in experiments. 

In the remaining of this section, we shall study the topological as well as 
analytic aspects of latitudinal states and spiral states. Analytic description 
of quasi-baseball states are more involved and will be discussed in Sec. [5] 

5.3 Hedgehog/Latitudinal States 

Using spherical coordinate system (9, (f>), we have a simple analytic form for 
latitudinal states 

'/'latitudinal = Ng dO , (75) 

where 9 is the polar angle, and Ng will be shown to be an integer below. 
The corresponding layer spacing can be calculated using Eq. ([12]) . with the 
metric tensor given by 

gee = R 2 , 94,4, = R 2 sin 2 6», g e4> = g^g = 0, (76) 

where R is the radius of sphere. We find that 

d = 2irR/\N \. 

Therefore Eq. (|75|) indeed describes a smectic state with equal layer spacing. 
The differential form Eq. (|75p is actually exact, with the corresponding phase 
field for the density wave given by: 

©latitudinal = Ng 9 + <9 . (77) 

Note, however, that '/'latitudinal is discontinuous at two poles, where 9 = 0, tt. 
The density profile nevertheless should be smooth everywhere. Using Eq. (fT]) 
we find 

dp(6,4>) = -pism.G{6)dG{0) 

= -piN e sin 0(6) d6. (78) 

which is nonsingular at both poles only if sin (9 = there. That is, 

6» = k tt, Ngir + Q = k' tt, (79) 

with k and k 1 are integers. Hence Ng must be an integer. In the setting of 
AB-block copolymer film, these boundary conditions imply that right at the 
center of disclinations, either density of specie^ or of specie B is maximum. 
The number of smectic layers between two poles is given by 

TT - / latitudinal = ^~Ng [ d9 = ^-Ng. (80) 

Zir J N l-K J 2 

Hence the geometric meaning of the integer charge Ng is twice the internal 
distance between two poles, i.e. the number of smectic layers between them. 
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This internal distance is independent of the contour connecting two poles, 
since ip is exact. Now consider an infinitesimal elastic deformation generated 
by a vector field X that preserves the locations of disclinations, i.e. X must 
vanish at both poles: 

X| w =X| s = 0. (81) 

The infinitesimal change of the differential form Sip induced by this elastic 
deformation is given by Eq. (|48]) : 



Therefore 



8ip = d$ = d(-ei x ip). (82) 



/ d(-eix.ip) = (-eixVOI 



s 

Sip 

N JN 

= -eX a iP a \ s + eX a iP a \ N =Q, (83) 



s 

N i 



where used was Eq. ([81]) . Therefore, 



S pS 

Sip = / Sip', (84) 



n Jn 



i.e., the internal distance Nqtt between two disclinations at poles is invariant 
under an elastic deformation that docs not displace two disclinations at poles. 
Ng is therefore indeed a topological quantity in the weak sense. 



5.4 Spiral States 

In the presence of two integer disclinations at both poles, a loop that encloses 
the north pole once is non-retractable, as long as the loop is not allowed to 
cut the disclination q The integral of ip along this loop therefore gives us 
an integer dislocation charge N^: 

ip = 2ttAV (85) 

For the latitudinal states, this integer is clearly zero, since ^latitudinal is exact. 
The geometric meaning of is the minimal number of layers intersecting 
a close loop winding around the north pole. A smcctic state with one +1 
disclinations at each pole, and with a dislocation charge associated with 
the loop 70, is a spiral state. Spiral states are therefore classified by an inte- 
ger dislocation charge N,;,. There are two equivalent ways to understand the 
integer N^. If we exclude both poles from our consideration, the resulting 
manifold is topologically a cylinder. It is then clear that Ng is the global 
dislocation charge associated with the non-retractable loop winding around 

22 Consider applying arbitrary elastic deformation that leaves disclinations fixed 
to such a loop. The resulting loop still encloses the north pole once. 
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the cylinder once. Equivalently, we may think of Eq. (|85|) as giving the dis- 
location charge of the defects at the north pole. Consequently, the defects at 
the north pole has both a disclination charge +1 and a dislocation charge 
A^, and therefore should be understood as a composite defects. Clearly the 
defects at the South pole is also a composite defects, with dislocation charge 
— N^. The total dislocation charge of a smectic orders on a compact substrate 
is always zero, according to the discussion following Eq. (|4Tj) . 

We can also obtain an analytic expression for spiral states in the region 
that is not too close to either pole: 

ip = N^d<j> + a e d6, (86) 



which clearly satisfies Eq. (|85|) . Note that ag can take arbitrary real value. 
The layer spacing for this state can be calculated using Eq. (fT2|) : 



27ri? 

d= . (87) 
y/a*+N*/ sin 2 9 

As long as |A^| -c \a$\ and sin 9 <C l-ZV© | / 1 a© | , this layer spacing is almost 
independent of 9, hence Eq. ([56]) is indeed a good approximation for the 
spiral states. On the other hand, when 6 — > 0, 7T, \tp\ — > oo, and the layer 
spacing calculated in Eq. ([57]) becomes very small. This means that Eq. {57J) 
is not a proper description of the spiral states near two poles. Indeed we 
will show below that the spiral states should really be understood as quasi- 
baseball states with two pairs of closely bound disclinations. To determines 
smectic layer configurations in regions sufficiently close to the poles, the 
exact locations of two +1/2 disclinations around each pole are needed. This 
information is however lost in the spiral state expression Eq. (|86|). leading 
to singular behavior of the layer spacing in Eq. (|87[). Finally when ag ~ N<f,, 
the description Eq. I|87p qualitatively breaks down everywhere. These states 
must be described as quasi-baseball states, which shall be studied in detail 
in Sec. [5] 

We have seen that in the case of latitudinal state Eq. ([75]) , the parameter 
Ng must be an integer, and is related to the internal distance between two 
poles. In the case of spiral state, however, this internal distance depends on 
the choice of contour and is not well defined. Therefore the corresponding 
parameter ag can take arbitrary real value. In fact, we can show that spi- 
ral states with the same A^ and different ag form a continuum, which are 
topologically identical to each other, and are mutually accessible by elastic 
deformations. To see this, let us consider a prototypical state 

V> = Nftfy, (88) 

which describes smectic layers in the direction of longitudes, shown in 
Fig. [T5K. Let us elastically twist the sphere around the z axis by a 9 depen- 
dent angle A<j) = — j^Q- A point with coordinates (6,<p) is deformed to be 

the point with coordinates {9',<j>') where 

9^9' = 9, 4> -> <{>' = <t> - ^-9. (89a) 
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Fig. 15 How to elastically cha nge the parameter ag for the spiral state Eq. (18611 , 
A) Starting from the state Eq. (|88[) with N$ = 2, we twi st the sphere by an total 
angle agir/Nj,. B) We end up with the spiral state Eq. (|86p . This process also shows 
that the internal distance between two poles in a spiral state is not a well defined 
quantity, and depends on the contour over which ip is being integrated. 

The new smectic state induced by this deformation, illustrated in Fig. [T5B. 
is then given by 



which is precisely the spiral state Eq. (|86p. expressed in terms of new co- 
ordinate (0',(j)'). Hence all spiral states with dislocation charge Ng can be 
continuously generated from the prototypical state Eq. i88\) by elastic defor- 
mations. By Eqs. (|89[) we also see that every time we twist the sphere by 
7r (that is, we swap the disclination pair in a counter clockwise direction), 
the parameter ag increases by N$. In the regime N$ <C ag, the layer spacing 
changes by a factor of N^/ag, with the resulting strain energy scaling as 



which only depends on but not on ag. Such twists of disclination pairs 
therefore costs low elastic free energy, as long as is of order of unity. 



6 Quasi-baseball States of Spherical Smectics 

Let us consider a quas-baseball state Fig. [TBR with four disclinations widely 
separated, and choose the coordinate system such that all four disclinations 
are on the equator. As illustrated in Fig. [17l we stereographically project 
the sphere from the south pole to the complex plane, such that the equa- 
tor is mapped into a unit circle. The resulting smectic pattern as well as 
disclinations a, b, c, d, ordered counter clockwise, are shown in Fig. I16B . Let 
us construct a loop (red) enclosing the ah segment of equator. The minimal 
number of smectic layers intersecting this red loop is N\. By deforming the 
loop to be infinitesimally close to the segment ab, one clearly see that N± is 
precisely twice the total number of smectic layers intersecting the segment 



*//(<?', <j>') = AV#(0', 0') = N+dtf + agd6', 



(90) 



E B ~BNM 



(91) 
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Fig. 16 (Color online) (A) A quasi-baseball state, with four +1/2 disclinations 
sitting on the equator. Each disclination core has one layer termination. Only two 
disclinations are visible here. (B) The same state stereographically projected onto 
the complex plane. The equator is mapped to the unit circle. The red contour 71 
encloses defects a, b, and intersects 12 smectic layers in total. The blue contour 72 
encloses defects b, c and intersects 8 layers in total. Hence this quasi-baseball state 
is labeled by two integer charges (Ni = 12, N2 = 8). 
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Fig. 17 (Side view) stereographic projection of a unit sphere onto a plane from 
the South pole. The equator is projected to the unit circle. 

ab of equator H Note also that because of the peculiarity of sphere, we can 
continuously deform the red loop, without cutting any defects, so that it 
tightly encloses the segment cd of equator. The integer N\ can not change 
during such a smooth process. Therefore Ni is also twice the total number 
of layers cut by the segment cd. Likewise, we can also introduce a (blue) 
loop enclosing the segment be of equator, which cut a minimal number N2 
of semctic layers. A similar argument then shows that N2 is also twice the 
total number of smectic layers intersecting the segment be, as well as by the 
segment da. We shall claim that two integers (Ni, N2) behave as global dis- 
location charges and completely determine the topological properties of a 
quasi-baseball state, under the assumption that disclinations arc immobile. 



23 In counting this number, a layer terminating/starting right at a or b is counted 
as half. 
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Let us set up the coordinate system in the complex plane such that the 
four defects are given by 

(o,6,c,d) = (l,e 4 * -l,e-*). (92) 

From the director pattern in Fig. [Til we easily deduce that all smectic layers 
intersect vertically with the equator. Let d be the equilibrium smectic layer 
spacing, which is uniform on the sphere, but might be slightly different with 
its energetically preferred value do, we find 

R(j) = N ± d, 
R{-k-4>)= N 2 d, 



From which we find 



Ni = <t> 

N 2 7T - (j) 

ttR 



Ni + N-2 



(93) 



The equilibrium position of defects are therefore completely determined by 
the two integers N± and N 2 . Large deviation of defects positions necessar- 
ily induces significant change of layer spacing, and therefore is energetically 
penalized. 

Consider a special case of quasi-baseball state with (j) = 0, and N\ = 0. 
The two +1/2 disclinations at z = 1 and z = e 1 ^ becomes degenerate and 
form one +1 disclination. This is clearly a latitudinal state with Nq = N 2 , 
now established as a special case of quasi-baseball states. Consider now the 
case 4> <^ 1, an d N\ <C N 2 . The two +1/2 disclinations at z = 1 and z = e 1 ^ 
form a close bound pair, the two at —1 and —e 1 ^ form clearly form another 
pair. If one realign the sphere such that two defects pairs are located at two 
poles, and calculate the integral of ip along a loop enclosing one defects pair, 
one finds the result 2nNi . This is a defining property of the spiral state with 
Ncf, = Ni, Eq. ([55]) . Hence spiral states are also special cases of quasi-baseball 
states. 

Eq. ([9"3"]) shows that N\ + N 2 of every low energy states is completely 
determined by the sphere radius R and the layer spacing d, while Ni can 
take arbitrary integer from to nR/d. Hence the total number of low energy 
smectic states on sphere is given by 

^(F fl <ifldg) = ^ + 1«^ > (94) 

i.e. scales as the square root of the total sphere area. 

Our purpose is to develop an analytic description ip for low energy smectic 
states that have no dislocation, i.e. dip = except at four disclination cores, 
and have equal layer spacing, i.e. 2it/\i(j\ = d « do- This turns out to be 
rather difficult. Let z = x + iy is the complex coordinates obtained from the 
stereographic projection. A complex differential form in the form of 



V/ = f(z)dz 



(95) 
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is called an analytic form, or holomorphic form, if f(z) is an analytic func- 
tion. Using the formalism we develop in Appendix [XJ one can show that the 
closed differential form ip that we seek for can be expressed as the linear 
superposition of the real part of an analytic form tJj c and an exact form da: 

ip(x,y) = Reip c (x,y) +da(x,y), (96) 
r(x,y) = f(z)dz, (97) 

where &{x, y) is a smooth (but not analytic) function on the xy plane. ip H = 
Reip c (x,y) is often called a harmonic form. Since an exact form da only 
describes an elastic deformation, the topological properties of the closed form 
if) is completely encoded in the analytic form f(z)dz. In particular, as we shall 
show in Appendix [AT the +1/2 disclinations of ij; are represented as poles 
and branch points of the analytic function f(z). We note, however, unlike ip, 
the layer spacing in the state described by ip H is generically not constant. 

We have therefore broken down our task of finding the analytic descrip- 
tions of quasi-bascball states into two steps: In the first step we find a har- 
monic form tp H that has the same topology as ip, but does not have equal 
layer spacing. In the second step we find the elastic deformation da that 
makes the linear combination ip H + da has equal layer spacing. In this work, 
we shall only focus on the first step, and leaving the second step to an sep- 
arate publication. There is a more physical way to understand this process: 
one can show that ip H describes the dislocation-free smectic ground state in 
the limit of vanishing layer compressional modulus B — > 0. Starting from this 
(rather fictitious) state and adiabatically tuning up the modulus B — > oo, we 
automatically obtain the real smectic ground state with nearly incompress- 
ible layers. 

The complex analytic form ip c = f(z)dz that describes a given distri- 
bution of disclinations is straightforward to find. As shown in Appendix [A| 
the winding number of f(z) at a point zq is just negative the disclination 
charge at the same place. Hence a form dz/ z describes a +1 disclination at 
the origin, while zdz describes a —1 disclination at the same point. A half 
integer disclination at z = a is naturally described by a multiple valued form 
dz/y/z — a. A quasi-baseball states with four +1/2 disclinations sitting at 
a,b,c,d respectively is described by the following 1-form with four branch 
points: 

Ae~' ia 

V-' c = . dz (98) 

(z — a)(z — b)(z — c)(z — d) 

Our task is then to determine the complex coefficient A e~ la in this holomor- 
phic form for a given quasi-bascball state with charges N± and A^. 

To avoid the ambiguity of sign associated with the multiple-valued func- 
tion in Eq. (|98p . we introduce two branch cuts along the equator, that con- 
nects branch points a and &, as well as c and d respectively, as illustrated in 
Fig. [TBI Ell The differential form Eq. (|98|) can be made single valued, as long 



24 There are topologically distinct ways to introduce branch cuts for a given pair- 
ing of branch points. Exactly which one we choose is however not important; what 
matters is that they serve as an impenetrable fence preventing us from tracing out 
a close loop that winds around a single branch point. 
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as it is not allowed to trespass the branch cuts. Furthermore, the resulting 
cut sphere becomes topologically indentical to a cylinder: there is only one 
nontractable loop, the solid loop 71 which encloses the branch cut ab coun- 
terclockwise. The integral of ip H along the close loop 71 defines an integer 
topological charge N±, which we introduced at the beginning of the current 
section: 



This integral is invariant under continuous deformation of the integration 
contour as long as it does not intersect either branch cut. By deforming the 
contour sufficiently close to the branch cut ab, one easily see that geomet- 
ric meaning of the integer charge N\ is twice the number of smectic layers 
intersecting the branch cut ab. 

We can also choose another equivalent branch cuts scheme, i.e. two equa- 
tor segments connecting ac and bd respectively, as shown in Fig. 1161 The 
resulting cut-sphere, topologically a cylinder, admits one nonretractable loop 
72, which encloses the branch cut ac counterclockwise, as well as encloses bd 
clockwise. The integral of ip c along 72 then defines the other integer topolog- 
ical charge N2, whose geometrical meaning is twice the number of smectic 
layers intersecting the branch cut ac (or bd): 



For given defects positions a, 6, c, d, two topological charges (N± , N2) uniquely 
determine, through Eqs. ([99]) and (jlOOp . the complex coefficient Ae~ la in the 
holomorphic form in Eq. (|98| . 

Another cquiavalent, but more elegant way of representing a multi- valued 
complex function on the Riemann sphere is to introduce a multi-sheets Rie- 
mann surface Each sheet is a copy Riemann sphere cut properly. Differ- 
ent sheets are inter-conncctcd appropriately in a way that is specified by the 
multi-valued function. 

The Riemann surface of the multi-valued function in Eq. (|98p consists 
of two spheres, one inside the other. As shown in Fig. [181 each of them is 
cut open through ab and cd respectively. E3 Finally two spheres are cross- 
connected along the branch cut ab and along cd respectively. The topology of 

25 Note that the same loop cut also encloses the other branch cut cd clockwise, 
since the infinity on the complex plane denotes the north pole on the sphere. 

26 We refer to Stone's notes [49| for an introduction of Riemann surface and mul- 
tiple valued function. 

27 Clearly we can also cut each sphere open along ac and bd, or ad and be. The 
topology of the final manifold obtained is however independent of this choice. 





(99) 




(100) 
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Fig. 18 The Riemann surface of the multiple valued function in Eq. (|98[) consists 
of two concentric spheres. The outer sphere is cut open in order to show the inner 
one. Both are cut along ab as well as along cd respectively. They are further cross- 
connected through the cuts. The resulting Riemann surface is topologically a torus, 
characterized by two non-retractable loops. 



Yi 




Fig. 19 Two ind ependent non-retractable loops 71 and 72 on the Riemann surface 
shown in Fig. [18] Solid lines are on the outer sheet, while dashed lines are on the 
inner sheet. The loop 71 lies on one single sheet and encloses branch points a and 
b. 72 lies on both sheets and encloses branch points b and c. Note that whenever a 
line intersect a branch cut, it leaves the sheet it resides and enters another sheet of 
the Riemann surface. 

the Riemann surface thus obtained is however a torus [49) , with two distinct 
nonretractable loops. On of them is 71 as shown in Fig. [19j which completely 
lies on one sheet of the Riemann surface. The other loop 72, also shown 
in Fig. [HI lies on both sheets of the Riemann surface. The integral of ip c 
along these two loops then define two dislocation charges through the same 
equations ((99)) and (|100|) . A curve that encloses one single branch point once 
is not a close curve in the Riemann surface. The curve encloses the branch 
point d (shown in the figure) twice is a loop that can be retracted to a single 
point. One can further show that the superposition of loops 71 and 72 can be 
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shown to be topologically identical to a loop enclosing branch points a and 
c only. 

The integral 

w = r 1 (z)= f - dt (101) 

J Z0 y/{t-a)(t-b)(t-c)(t-d) 

is usually called a an elliptic integral. It defines w as a multi- valued function 
of z, because the value of the integral depends on the path connecting zq 
and z. Every time the path wraps around the non-retractable loop 7, once, 
w changes by an amount Wj, given by 

dt i = 1,2. (102) 

y/(t-a)(t-b)(t-c)(t-d)' 

w = I^ 1 (z) is therefore only determined by z up to integer multiples of ui\ and 
L02- The inverse function z = I(w) is nevertheless single- valued and is called 
an elliptic function, which has two periods uo\ and u>% given by Eq. (|102|) : 

z(w + lui) = z(w + L02) = z(w). (103) 

Remember that the physical meaning of ip = Re tp c is the differential of 
the local phase field of the smectic density modulation, therefore the integral 
of Re ip c is the multiple- valued phase field 0{z), which, according to Eqs. (|99l 
I lOOp . may have uncertainty 27riVi or 2irN2- Multiplying a constant Ae~ la to 
Eq. (jlOip and taking the real part we find 

0(z)=Ref V c = Re (Ae _<a J _1 («)) . (104) 

J Zq 

Comparing Eq. fTU2]) with Eq. ||j55J) and pU0|) . we find that the charges Ni 
of the smectic state ip c are related to the complex periods Ui of the elliptic 
integral via 

2tt N t = Re (tOi A e~ ta ) , i=l,2. (105) 

This establishes the relation between the phase field &(z) for spherical smec- 
tic and the elliptic integral Eq. (|101[) . 

For four disclinations on equator as given by Eq. (|92|) and for given two 
charges (iVi,.^), let us solve Eq. (f99]) and (| 1 00[) for two parameters A and 
a. Parameterizing z = exp iip along the integrating contour, the two integrals 
Eqs. (|99H100p can be transformed into 

. i> c = iAe-^- la dif (106) 

71 Jo v sin ip sm{<f> — ip) 

= Ae -i4>-ia I dL P — ( 107 ) 

7 2 J 4, i/sin< / 9sin((/3 — 0) 
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Fig. 20 Starting from a state (Ni = 2, N2), we twist the sphere counter-clockwise 
by 7r so as to exchange the defects pair a, b, we obtain a state (TVi = 2, iV2 + 
Ni). These two states are therefore topologically equivalent. Note that this twist 
operation changes the layer spacing by a factor of N1/N2. 



These two integrals can be expressed in terms of complete elliptic integrals 
of the first kind: 



2ttN 1 = 2A K \ | J sin + aj , (108a) 
2ttN 2 = 2AK ( cos| ) cos ( | + a] . (108b) 



Derivation of these results is relegated to Appendix [Bl Solving Eqs. (|108|) for 
A and a we find 

Ae^=(-^-i**)e«<>, (109) 
^A-(cos|) K (sin I) ) 

where A'(-) is the complete elliptic integral of the first kind. 



6.1 Purely Topological Classification 

Spiral state are especially interesting because of their closely bound defects 
pairs. As illustrated in Fig. [20l let us start with a state with charges (iVi = 
2, N2), fix the south pole and elastically twist the whole sphere around the z 
axis by an angle of 7r so that the defects pair a, b at the north pole exchange 
their positions. Note that all smectic layers also deform with the twist Fi . 
The branch cut which connects a and d before the twist operation, shown 
blue in Fig.[20l\. after the 1 twist operation, connects b and d intead, shown red 
in Fig.[20t3. The corresponding integer charge associated with this branch cut 
is N2 before the twist, and becomes N2 + Ni after the twist. We therefore 
conclude that counter-clockwise twist of the defects pair a, b leads to the 
following transformation of the dislocation charges 

{N 1 ,N 2 )^(N 1 ,N 2 + N 1 ). (110) 

28 It is amusing to point out the similarity between this twist and a position 
exchange of an anyon pair. 
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A B 

Fig. 21 A): A unit cell in the comp lex p lane with bases periods u>i and W2 for the 
elliptic function I(z) defined in Eq. (|101[) . can be glued into a torus. B): Similar to 
a crystal, we have freedom in choosing the two bases periods to be (lji,^), or to 
be (lui, uii + lu 2 )- 

Obviously clockwise twist of the pair (a, b) leads to the transformation of 
dislocation charges: 

(iVi,iV 2 ) ->■ (iVi,7V 2 -iVi). (Ill) 

We have then proved that, as long as defects are allowed to move, states 
(Ni, N 2 ± Ni) are topologically identical to the state (Ni, N 2 ). 

More generally, we may try to twist the defects pair by an arbitrary angle 
6, so that the defects configuration in the twisted state does not exactly 
coincide with that of the initial state. As long as we do not focus on the 
smectic layer configuration near the polar regions, however, we are not able 
to defects the change of defects configuration. Equivalcntly we are not able 
to unambiguously choose the branch cut that connect one of the defects near 
the north pole to one near the south pole. This leads to a complete loss of the 
quantization of the charge N 2 — it can take arbitrary real value. We therefore 
end up with a more coarse-grained description of spiral states as discussed 
in Sec. Ed e.g. Eq. ®. 

Note, however, our twist manipulation of defects is topologically applica- 
ble for arbitrary states, with Ni or N 2 not necessarily small. Starting with an 
arbitrary state with N\ > iV 2 , and twist an appropriate defects pair so as to 
obtain a pair (N[ = N x - N 2 , iV 2 = N 2 ). If N[ > 7V 2 still holds, we simply re- 
peat the same twist operation. If, after the twist, we find N 2 > N[, we twist a 
different pair so that we have a new pair (N[, N^ — N[). This process must end 
after finite steps, where we end up with one charge vanishing. Note that this 
reduction process is precisely the Euclidean algorithm for finding the greatest 
common divisor gcd{N\, N 2 ) between two integers (Ni, N 2 ). Therefore at the 
final state we find must have charges (gcd(N\, N 2 ), 0), which is a latitudinal 
state. Hence we have shown that if disclinations are allowed to move, all 
quasibaseball states are topologically identical to some latitudinal state, with 
only one nonzero integer charge. This settles the problem of classifying all 
dislocation-free smectic states on sphere based on purely topological consid- 
eration. We note however, this classification has little physical relevance, due 
to the significant change of layer spacing caused by twist processes. 

We conclude this section with an alternative interpretation of the trans- 
formations Eqs. (|1 lOtil 1 1|) . We have already seen that Eq. ()101|) defines a 
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doubly periodic function on the complex w plane. This doubly periodicity 
defines a lattice structure in the w complex plane with a unit cell spanned by 
two complex numbers lo\ and u>2, as shown in Fig. I21I A We have freedom in 
choosing the two periods of a doubly periodic function. As shown in Fig. l21B . 
if (wi,u>2) are two periods and define a unit cell, then so are (wi,u;2 + u>i)- 
In fact for arbitrary integers k, I, m, n with kn — Im = ±1, 

lu[ = ku>i + I0J2, 
uj' 2 = mu>i + nw2, 

are also periods and form a unit cell. This is because the u)\ and u>2 can 
also be expressed as integer combination of u>[ and lo' 2 . This redefinition of 
unit cell is of course well known in study of crystal structures. Twisting of a 
defects pair by n therefore is equivalent to a redefinition of the unit cell in 
the w plane for the doubly periodic function z = I(w). 

7 Conclusion 

Smectically ordered systems confined on a compact curved substrates, such 
as torus and sphere, exhibit rich topological as well as geometric structures. 
The formalism of differential forms provides a natural framework for analyses 
of these structures. 
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A Complex Differential Forms and Holomorphic Forms 

Our goal is to develop an analytic description in terms of differential forms for 
smectic states with four +1/2 disclinations on a sphere. The major difficulty of this 
problem lies in the fact that, in the presence of +1/2 disclinations, it is impossible to 
choose a sign for the differential form ip without any ambiguity. Parallel transport 
of ip (or the director field n) around one half integer disclination yields —ip (or 
— h). Hence we expect that smectic orders in quasi-baseball states are described by 
multi-valued differential forms. Our formalism can be applied to arbitrary curved 
surface, instead of sphere only. A basic introduction o f re levant formalisms can be 
found in chapter 7, Vol II of M. Stone's lecture notes [49] . 

In this section we shall exclusively use conformal coordinate system, where the 
metric tensor is isotropic, i.e. 

9ab = ^gS ab , g ab =-^=S a b- (112) 
v 9 

Note that g is generically a function of coordinates. Gauss proved that conformal 
coordinate system always exists for arbitrary two dimensional surface. Let 

ip = p(x,y)dx + q(x,y)dy (113) 

be a differential form expressed in conformal coordinates (x,y), the Hodge dual of 
the 1-form ip, denoted as -kip, is defined to be 



*lj} = -q(x, y)dx + p(x, y)dy. 



(114) 
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A more complete definition of Hodge dual using arbitrary coordinate system is 
given in Appendix [C] A harmonic form is closed form whose Hodge dual is also 
closed, i.e. 

dip = d*ip = 0. (115) 

It is easy see that 

dip = (d x q — d y p)dx A dy, (116a) 
d-kip = (d x p + d y q)dx A dy, (116b) 

Therefore ip is harmonic if and only if 

d y p = d x q, d x p = -d y q, (117) 

which are precisely the Cauchy-Riemann c ondit ions for a complex analytic function 
f(z) — p(x,y) — iq(x,y). Therefore Eq. (|113|l is a harmonic form if and only if 
p(x, y) and —q{x, y) form a conjugate pair of harmonic functions. This explains the 
terminology "harmonic form". 

If Eq. (|113[1 is harmonic, the complex function f(z) = p{x,y) — iq(x,y) is 
analytic, i.e. holomorphic, where z = x + iy. We can therefore introduce an analytic 
1-form 

ip c = f(z)dz = (p-iq)dz, (118) 
which depends only on z, but not on z. Defining the complex conjugate 

i> c — (p + iq)dz, 

we find the real harmonic form ip and the complex holomorphic form tf) c are related 
to each other via 

i> = \{^ C + -¥) =Rei/j°, (H9a) 
ip c = ip + i * ip. (119b) 
The dislocation charge associated with a loop 7 is then given by 

ip = Red) ip c . (120) 

J-y 

Using the conformal metric, the magnitude of ip can be calculated easily: 

\ip\ 2 =g ah ip a ip b = ^-{p 2 + q 2 ), (121) 

which is related to layer spacing by Eq. (|12|l . The associated normalized 1-form 
(with unit magnitude) therefore equals to the nematic director form: 

h = \ip\~ ip = V (pdx + qdy) 

vp + <r 

= h x dx + hydy. (122) 



Let us define the magnitude of the corresponding complex form ip c Eq. (|118 



|^| 2 = ^=|/(*)f = ^(P + 1) = l^l 2 , (123) 
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hence is the same as the magnitude of the real form if). The normalized complex 
1-form is then given by 

h c = -^—ih c = , V ^ (p - ig)dz (124) 



Comparing Eq. (|122p wi th Eq , (JT24JI, we find that, if the smectic state is described 
by a harmonic form Eq. (|113[) . the phase of complex analytic function f(z) = p — iq 
encodes the direction of the local nematic order relative to the conformal coordinate 
system (x,y), that is 



Arg/O) = -tan" 1 (125) 



Let us consider the nematic director Eq. (|122[) as a function of coordinate z. If 
there is a nematic disclination at point zo, the nematic director winds around the 
unit circl e once as z winds along a loop enclosing the point zo couter-clockwise. 
Eq. (|125[) then tells us that the phase of the analytic function f(z) winds the 
complex unit circle once clockwise. More generally, the winding number of the 
analytic function f(z) along a close loop is negative the disclination charge enclosed 
by this loop. In other words, positive (negative) nematic disclinations are described 
by poles (zeros) of the meromorphic function f(z). This relation between phase of 
analytic function and two dimensiional orientational orders has been known for a 
long time and pro vides t he basis for the application of complex analysis in nematic 
orders on sphere |29H53j . We also present a simple derivaion of this relation using 
differential forms in Appendix [D] The basic method works, however, on arbitrary 
curved surface, as long as a conformal coordinate system is used. 

Our purpose is however to use analytic (excluding possible singularities) func- 
tion to describe smectic ordering on curved surface with no dislocation except at 
the disclination cores. For this purpose we need to establish the relevance of har- 
monic forms, which is equivalent to holomorphic forms, as we have shown above, 
to smectic ordering on curved substrates. This is achieved by the Hodge decompo- 
sition theorem of differential forms 0: On a compact manifold, an arbitrary closed 
1-form ip can be decomposed into an exact form da, and a harmonic form tp H , 

ip = da + ip H . (126) 

Remembering that an exact form da describes an elastic deformation and does not 
change any topological property. Let ip describes the smectic state that we aim to 
study, the Hodge's decomposition theorem says that its topological properties can 
always be captured by a harmonic form tp . In particular, the integrals of ip and 
tp H along arbitrary close loop equals to each other: 

i II 



V> = ^V> , (127) 

since the integral of an exact form da along a closed loop 7 always vanishes. Hence 
we conclude that in order to understand the topological structures of smectic or- 
dering on curves surface, we only need to study harmonic forms on sphere. 



A.l Application to Sphere: Hedgehog States 

Let us apply the formalism of holomorphic forms to sphere. We first stereogra phi- 
cally project a unit sphere into a plane from the south pole. As illustrated in Fig. 1171 



What we have invoked here is actually only a special case of the Hodge decom- 
position theorem for arbitrary differential forms. 
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an arbitrary point (6, cj>) on the sphere is projected to the point (r = tan#/2,0) 
using polar coordinate system. The north pole is projected to the origin, the south 
pole to the infinity, while the equator 6 = tt/2 to the unit circle r = 1. Introducing 
the Cartesian coordinate (x, y) in the plane: 

x — r cos0, y — r sin</>, (128) 

we easily verify that they are indeed conformal coordinates: 

9xx = 9yy = M + r 2)2 ' Sxy = 0. (129) 

We can then subsequently form a complex coordinate z = x + iy. 
Let us first consider the holomorphic 1-form 

4> c = f(z)dz = -dz. (130) 

z 

The winding number of - around Zo — is clearly — 1, hence Eq. (|130p describes a 
+1 nematic disclination at z = 0, i.e. the north pole, according to the stereographic 
projection. In order to look at the configuration around the south pole, we should 
project from the south pole. The resulting new complex coordinate (also conformal) 
w is related to z via w = 1/z. In terms of w the 1-form tp c is given by 

Tp c — — dz = — — dw. (131) 

z w 

Hence the winding number around w — (z — oo) is also —1, which corresponds 
to a +1 nematic disclination at south pole. This is of course rather expected, as 
the total disclination charge has to be +2 by the Gauss-Bonet-Poincar e theorem. 

Let us look at the corr espon ding real 1-form, using polar expression z = re"* 
as well as the relation Eq. (|119|l we find 

ip = Reip c = -dr = d (log r) = d {log tan Q , (132) 

which is actually an e xact form. The layer spacing can also be calculated, using 
Eq. (H2J) and Eq. lfl2"3]) , as well as Eq. Ip9]) . 

j 27r 2nr 

W\ = T+r 5 ' (33) 

which approaches zero as r — > (north pole) and r — > oo (South pole). Hence 
smectic layers become infinitely dense around both poles. This peculiarity is a 
consequence of -0 being harmonic. As it stands the harmonic form tf> does not 
describe the ground state with minimal energy. The groun d sta te which clearly 
prefers equal layer spacing, should be solved by using Eq. (|126(l and minimizing 
the total elastic free energy over the elastic deformation da. Such an analysis is left 
to a separate publication. 

Let us integrate ip and find the phase field 0: 

tan | , % 

0(6) - O = log — £ (134) 
tan 

where 80 and 6>o are two arbitrary constants. Hence all smectic layers form circles 
of co nsta nt latitude, while nematic director is parallel to lines of longitude. Hence 
Eq. (|130[) describes a latitudinal state. This result can also be seen by calculating 
the dislocation charge associated with the loop 7^: 

ij) = Re<f ip c = 0, (135) 

'14, J 14, 

which corresponds to a latitudinal state. 

Finally we note that description of smectic layers in terms of 1-form Eq. (|130[) 
only works outside the disclination core regions. 
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A. 2 Spiral States and Composite Defects 

Let us now consider a slightly more complicated 1-form: 

ip c = f(z)dz = A e~ ia -dz. (136) 
z 

It still describes two +1 disclinations sitting on two poles. On the other hand, the 
real 1-form is given by 

ip = Re(ip c ) = A cos a —dr + A sin a d<f>, (137) 
r 

where </> is the polar angle of the spherical coordinate system. Integrating ijj along 
a loop 70 enclosing north pole once, one must obtain an integer multiple of 2n, if 
layers are well defined. That is, 

ip = 2n A sin a = 2fc7r, (138) 



which imposes constraint on the parameters A and a. Eq. (|136|) therefore describes 
a spiral state with global dislocation charge (integer) given by A sin a. The discli- 
nations at both poles therefore are composite defects with both disclination charge 
and dislo cation charge. It is rather interesting to see t hat th ese composite defects 
Eq. (|136|l can be generated from pure disclination Eq. (|130|l by simple multiplica- 
tion of a phase factor. 

It is also easy to find the locus of constant phase cj>o of density wave by inte- 
gration: 

log r + cf> tana = (139) 

which is a logarithmic spiral. Again we note that description in terms of differential 
form only works outside the disclination cores, i.e. tan -1 ao < r < tan -1 I/do- 

More generally, let ijj c = f(z)dz be a holomorphic (except a few poles or branch 
cuts) form describing some particular smectic state with no dislocation. The integral 
of ip c along a close contour 7 is generically a complex number: 



tP c = * f(z)dz = Br+iBi. (140) 
7 Jo 

Using Eq. lfTT9) l we find that 

ip = B R , <b *ip = Bi. (141) 



That is, Br and Bi are the dislocation charges assoc iated with the loop 7 for 
smectic states ijj and *V respectively. We note that Eq. (|140p is the sum of residue 
of ip = f(z)dz enlcosed by the loop. 

Let us apply a global rotation of the form i{j c by multiplying a complex phase 
factor: 

i> c -> i> c ' = e- ia ip c . (142) 
The new real differential from ip' is then given by 



ip 1 = Re ip c = cos ai/i + sin a -kip, 



(143) 
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where -kip is the Hodge dual of ip, defined by Eq. (|114p . The dislocation charge 
associated with tp' is then given by 

f tp' = Re (f> tp c = B R cos a + Bi sin a, (144) 
Jc J c 

which is generically different from Br. Note that \ip c \ = \tp c '\, which means that 
layer spacing are the same in two states tp c and tp°'. Multiplication of a complex 
phase factor to a holomorphic form therefore provides a basic mechanism for c hang- 
ing dislocation charge without changing layer spacing. Note howeve Eq. (|144[) must 
be integer multiple of 2-n. This means that for given value of Br and Bi, only 
discrete choices of the angle a are compatible with smectic ordering. 



B Calculation of Integrals Eqs. ()106H107p 



To calculate the integral in Eq. (|106[1 . we use the identity 

1 ,x 1 



sm 



(ip) sin(0 — (p) = — cos(2(/3 — <j>) cos <f> (145) 



and introduce a new integral variable j3 = 2ip — tp. The original integral is then 
transformed to 

tp — — — e 



\[2 J -4, Vcos (5 — cos p 

iA-y-i* r* dp (146) 

2 J-<t> . /air, 2 t - a ;r,2 £ 



r<t> 

Introducing a new variable t through 



sm 0/2 



we transform the integral Eq. (|146|) into a complete elliptic integral of the first 
kind: 

4> c = 2iAe-^- 



, f 1 d 



)(l-t 2 sin 2 |) 



2iAe'i' t '~ la K (sin^ ) . (148) 



2, 

To calculate integral Eq. (|107p . we make variable transformations: 

<fi — TV — Lp, (j> — TV — 0. (149) 

Eq. (|107|l then reduces to 

ip c = Ae-^- ia [ d ^ (150) 



\/sm(psm(4>— tp) 

which as the same functional form as Eq. (|106p . Now using the result Eq. (|148|l we 
find 

C C = 2Ae-^- la K fsin^ 

V 2 

= 2Ae-^- ia K (cos I ) . (151) 
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Combining Eqs. (fTIHj) and (fT5T|) with Eqs. (|99lll00p we finally obtained two 
equations that uniquely determine A and a: 



2-kN 1 = 2AK ( sin | 1 sin ( | +a ) , (152) 



2wN 2 = 2A K ( cos |1 cos I | + a ) . (153) 



C Hodge Dual and Codifferential 

In this appendix, we list the relevant results about Hodge dual of differential forms 
as well as codifferential operator. All these resu lts can be found in standard text- 
books of differential geometry, see, for example [35U13II49] . 

In a d dimensional manifold, the Hodge dual of a p-form is a (d— p)-form defined 

as 

★ dx 11 A ■ ■ ■ A dx lp = n 1 yfgg njl ■ ■ ■ g lpJp (154) 



(d-p)\ 

C . . flrp3p + l . . . Arf^d 

lz 31'--JpJp+l'--Jd u ' J ' U " L ■ 

For a two dimensional Riemannian manifold with orthogonal coordinate system 
(x, y), we have 

★ 1 = ^fgdx A dy, (155a) 
-kdx = yfgg xx dy, (155b) 
*dy = — yfg g vv dx, (155c) 

■kdx Ady = (155d) 
V9 

It is clear that when acting on a p-form twice, we have 

*W = (-l) p(d - p) uj p . (156) 
For a conformal coordinate system, g xx = g vv — 1/y/g, and Eqs. (|155|l reduces to 

★ 1 = yfgdx A dy, (157a) 
■kdx = dy, (157b) 
*dy = —dx, (157c) 

■kdx Ady = — . (157d) 
V9 

If a and /? are two p-form fields, a A-kb is a d-form field, which can be integrated 
over the whole manifold. This allows us to define an inner-product of two p-form 
fields 

(a,/3) = {/3, a) = f aA*f3 



■i J dwa 4l ...< p /9 il - < » ) (158) 



where in the last equality, the indices of j3 have been raised using the metric tensor. 
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The codifferential operator 5 is defined as 

5 = (-l) dp+d+1 *d* (159) 
— > — * d-k, for d = 2. 

Note that 5 maps ap form into a (p— l)-form, i.e. it works in the opposite direction 
as d. 5 acting on 0-form is defined to be zero. It can be shown that 5 is the adjoint 
of d, i.e. 

{a,dp) = {5a,P). (160) 

One can also show that 

5 2 = d 2 = 0. (161) 
The self-adjoint Laplace-Beltrami operator acting on differential form is defined 

as 

Aa = (Sd + dS)a. (162) 

One can show that 

(a, A3) = (-4a, /3) = (da, dp) + {5a, 5/3). 

(163) 

Setting a = (5, one find 

(a, Aa) = (da, da) + (5a, 5a) > 0. (164) 
Hence the Laplace-Beltrami operator must have non-negative eigenvalues. 



D Differential Form and Nematics on Curved Substrates 

In this appendix we present a simple by elegant derivation of the connection be- 
tween equilibrium nematic texture on curved substrates and holomorphic functions 
using differential forms. Let D a denote the covariant derivative associated with the 
symmetric Levi-Civita connection of the manifold, we easily see 

D a hp — Dpn a = dahp — dph a (165) 

Consider the differential form h = h a dx a describing the nemati c dir ector field, and 
for simplicity using a conformal coordinate system. Using Eq. (|158[) we find 

dh = (difi2 — d2"hi)dx 1 A dx 2 , (166) 
(dh,dh) = J dcL>g^ 1 (dih 2 — c^m) 2 . (167) 

We recognize that this is precisely the in-plane bending component of the Frank 
free energy for nematic director in curved space. On the other hand, we also have 

5n — -kd-k (hidx 1 + fi2dx 2 ) 
= -kd(fi\dx 2 — ^dx 1 ) 
= * (d\h\ + 82112) dx 1 A dx 2 

= -^p (dini + d 2 h 2 ) (168) 
y9 
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One can then show that 

(8h,8n) = j du{D a n° ! ) 2 , (169) 

is the splay component of the Frank free energy. 

The intrinsic part of the Frank free energy for nematic order on curved 
substrate then can be expressed in terms of differential form as 

F A = ^-(8n, 8ti) + ^j-{dn, dn). (170) 

When K\ = K3, this can be further simplified using the Laplace-Beltrami operator: 

F A = ^(n,Ah). (171) 

Minimizing the Frank free energy over n subject to the constraint \n\ = 1, one 
find the Euler-Lagrange equation for n: 

Ah + X(x)h = (Sd + d5)h + X(x)h = 0, (172) 

where X(x) is a Lagrange multiplier field which generically depends on x. Let the 
function (0-form) <f>(x) satisfy the same equation as h, i.e. 

A(p + X<p = 6d(f>(x) + X(f>(x) = 0. (173) 

we can show that the 1-form tp = (ph is a harmonic form, i.e. 

Aip = {8d + dS)tp = 0. (174) 

Hence we have 

(p{x) 

Conversely, if ip(x) is a harmonic form, and let us define h(x) = ip(x) /\ip(x)\ 
(hence h(x) is properly normalized), and X(x) = —8d\ip(x)\/\ip(x)\, then one can 
show that h(x) automatically satisfies Eq. (|172[) . i.e. it indeed describes an equilib- 
rium nematic texture. 

This establishes the connection between the equilibrium nematic texture and 
harmonic form tp(x). Using the relation between harmonic form and holomorphic 
form, we further obtain the (well known) relation between equilibrium nematic 
texture and holomorphic functions. 
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